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A linear,  inviscld  theory,  termed  the  coupled  basins  theory,  has  been 
developed  to  analyze  the  response  to  periodic  incident  waves  of  an  arbitrary  shape 
harbor  containing  several  interconnected  basins.  The  region  of  consideration  is 
divided  into  an  open-sea  region  and  several  inner-basin  regions  (the  number  depending 
on  the  harbor  geometry).  The  solution  in  each  region  is  formulated  as  an  integral 
equation  in  terms  of  the  normal  velocity  at  the  entrance  and/or  at  the  common  boun- 
daries between  regions.  An  approximate  method  is  used  to  solve  the  integral  equation 
by  converting  it  to  a matrix  equation.  The  initially  unknown  boundary  condition  at 
the  entrance  is  determined  by  matching  the  wave  amplitudes  and  their  normal  deriva- 
tives at  the  harbor  entrance  and  at  all  the  common  boundaries.  The  solution  for  the 
response  and  the  amplitude  distribution  within  the  complete  harbor  can  then  be  ob- 
tained. It  has  been  found  that  the  coupled-basins  theory  gives  results  which  agree 
well  with  experiments  both  for  an  irregular  shape  harbor  as  well  as  for  a harbor 
composed  of  two  connected  circular  basins.  Various  aspects  of  the  response  of  harbors 
composed  of  several  types  of  circular  connected  basins  as  well  as  circular  harbors 
with  rectangular  entrance  channels  have  been  investigated.  It  is  found  that  to  a 
first  approximation  the  response  of  a coupled  harbor  system  can  be  constructed  by 
superposing  the  response  of  the  individual  harbors.  Certain  aspects  of  the  effect 
of  viscous  dissipation  on  harbor  resonance  are  discussed.  Some  attention  is  given 
to  problems  of  scaling  model  results  to  the  prototype  harbor. 
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ABSTRACT 


A linear,  inviscid  theory,  termed  the  coupled  basins 
theory,  has  been  developed  to  analyze  the  response  to  periodic  incident 
waves  of  an  arbitrary  shape  harbor  containing  several  interconnected 
basins.  The  region  of  consideration  is  divided  into  an  open- sea  region 
and  several  inner -basin  regions  (the  number  depending  on  the  harbor 
geometry).  The  solution  in  each  region  is  formulated  as  an  integral 
equation  in  terms  of  the  normal  velocity  at  the  entrance  and/or  at  the 
common  boundaries  between  regions.  An  approximate  method  is  used 
to  solve  the  integral  equation  by  converting  it  to  a matrix  equation.  The 
initially  unknown  boundary  condition  at  the  entrance  is  determined  by 
matching  the  wave  amplitudes  and  their  normal  derivatives  at  the 
harbor  entrance  and  at  all  the  common  boundaries.  The  solution  for 
the  response  and  the  amplitude  distribution  within  the  complete  harbor 
can  then  be  obtained. 

It  has  been  found  that  the  coupled -basins  theory  gives  results 
which  agree  well  with  experiments  both  for  an  irregular  shape  harbor 
as  well  as  for  a harbor  composed  of  two  connected  circular  basins. 
Various  aspects  of  the  response  of  harbors  composed  of  several  types 
of  circular  connected  basins  as  well  as  circular  harbors  with  rectan- 
gular entrance  channels  have  been  investigated.  It  is  found  that  to  a 
first  approximation  the  response  of  a coupled  harbor  system  can  be 
constructed  by  superposing  the  response  of  the  individual  harbors. 
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Certain  aspects  of  the  effect  of  viscous  dissipation  on  harbor 


resonance  are  discussed.  Some  attention  is  given  to  problems  of 


scaling  model  results  to  the  prototype  harbor 
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CHAPTER  1 
INTRODUCTION 

Waves  which  propagate  from  the  open-sea  into  bays  and  harbors 
usually  experience  changes  in  wave  direction,  shape,  and  amplitude 
due  to  the  local  bathymetry  and  horizonal  geometry  of  the  embayment. 

The  harbor  shape  and  the  internal  reflections  associated  with  the  harbor 
boundary  can  cause  amplification  or  attenuation  of  the  incident  wave 
system.  This  phenomenon  usually  is  referred  to  as  harbor  resonance, 
seiche,  or  harbor  surging  and  can  be  the  cause  of  significant  damage  to 
moored  ships  and  adjacent  structures  especially  if  the  resonant  period 
of  the  ship-mooring  system  is  close  to  that  of  a mode  of  oscillation  of  the 
harbor.  In  addition,  the  currents  induced  by  these  harbor  oscillations 
can  create  navigational  hazards  near  the  entrance  and  within  the  harbor. 

For  an  existing  harbor  that  experiences  operational  problems  due 
to  the  effect  of  resonant  oscillations,  corrective  action  must  be  taken  to 
reduce  or  eliminate  such  oscillations.  In  the  case  of  new  construction  an 
attempt  should  be  made  to  use  a harbor  geometry  that  will  be  free  from 
possible  resonance  effects  for  the  local  wave  conditions.  In  either  case 
an  analytical  method  for  predicting  resonance  and  the  associated  wave 
amplitudes  is  quite  useful  for  preliminary  investigations  as  well  as 
providing  a guide  for  experimental  studies  if  they  are  deemed  necessary. 
The  existing  methods,  to  be  discussed,  are  not  advanced  to  the  state 
where  the  resonance  characteristics  of  a harbor  can  be  determined 
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accurately  analytically;  however,  as  mentioned,  these  methods  certainly 
can  assist  the  engineer  in  early  design  stages,  in  assessing  existing 
problems,  or  in  conducting  experiments  in  the  laboratory. 


1.  1 PREVIOUS  HARBOR  RESONANCE  STUDIES 

Previous  analytical  and  experimental  studies  of  harbor  resonance 
can  be  divided  into  two  groups.  The  first  deals  with  the  problem  of  wave- 
induced  oscillations  in  harbors  of  simple  geometry  such  as  circular, 
rectangular,  or  combinations  of  these  simple  shapes,  e.  g.  , McNown 
(1952),  LeMehaute  (1961),  Miles  and  Munk  (1961),  Ippen  and  Goda  (1963), 
Raichlen  and  Ippen  (1965),  Miles  (1970),  and  Carrier,  Shaw,  and  Miyata 
(1971).  The  second  group  of  investigations  is  concerned  primarily  with 
harbors  of  complicated  geometry,  e.  g.  , Knapp  and  Vanoni  (1945), 

Wilson  (1959,  1960),  Wilson,  Hendrickson  and  Kilmer  (1965),  Leendertse 
(1967),  Hwang  and  Tuck  (1970),  Lee  (1969),  and  Lee  and  Raichlen  (1970). 
In  this  section  the  work  of  these  investigators  will  be  briefly  discussed 
in  the  order  just  presented. 

McNown  (1952)  investigated  some  of  the  response  characteristics 
of  a circular  harbor  of  constant  depth  excited  by  wave  incident  upon  a 
small  entrance  gap.  In  his  analysis  it  was  assumed  that  the  crest  of  a 
standing  wave  (antinode)  occurred  at  the  entrance  when  the  harbor  was 
in  resonance.  A similar  method  was  applied  to  rectangular  harbors  by 
Kravtchenko  and  McNown  (1955).  Thus,  for  resonant  motions,  this 
assumption  led  to  a boundary  condition  identical  to  that  for  a completely 
closed  basin.  Therefore,  the  wave  periods  associated  with  resonant 
oscillations  would  be  those  which  correspond  to  the  eigenvalues  for  the 
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complicated  shape.  In  the  case  of  a rectangular  harbor,  Miles  and  Munk 
applied  this  general  formulation  and  found  for  periodic  waves  that 
narrowing  the  harbor  entrance  leads  to  an  enhancement  of  harbor 
surging  (oscillation)  instead  of  a reduction.  This  phenomenon  was 
termed  the  "harbor  paradox".  The  increasing  wave  amplification  at 
resonance  as  the  entrance  width  decreased  was  also  found  by  LeMehaute 
(1955)  and  Ippen  and  Goda  (1963). 

Ippen  and  Goda  (1963)  employed  the  Fourier  transformation  method 
to  evaluate  the  waves  radiated  from  the  entrance  of  a rectangular  harbor 
to  the  open-sea  and  the  method  of  separation  of  variables  for  the  region 
inside  of  the  rectangular  harbor.  A solution  for  the  response  was  obtained 
by  matching  average  amplitudes  from  the  two  regions  at  the  entrance. 

They  found  good  agreement  between  experiments  and  the  theory. 

Ippen  and  Raichlen  (1962)  and  Raichlen  and  Ippen  (1965)  have  studied, 
both  theoretically  and  experimentally,  the  wave -induced  oscillations  in  a 
rectangular  harbor  connected  to  a larger  highly  reflective  rectangular  wave 
basin.  They  found  that  the  response  characteristics  of  the  harbor  were 
radically  different  from  a similar  prototype  harbor  connected  to  the  open- 
sea.  The  response  curve  of  the  former  was  characterized  by  a large 
number  of  closely  spaced  resonant  maxima  where  for  the  harbor  connected 
to  the  open-sea  only  several  distinct  resonant  modes  of  oscillation  would 
be  experienced  for  the  same  wave  period  range.  It  was  pointed  out  that 
to  reduce  the  coupling  effect  of  the  reflection  of  the  wave  energy  which  is 
radiated  from  the  harbor  entrance  efficient  wave  absorbers  were  necessary 
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in  the  main  wave  basin.  A subsequent  study  by  Ippen,  Raichlen,  and 
Sullivan  (1962)  showed  that  the  coupling  effect  is  indeed  significantly 
reduced  by  wave  absorbers  in  the  main  wave  basin. 

Most  recently.  Miles  (1970)  re-examined  the  "harbor  paradox" 
using  electrical  circuit  analogy  and  used  as  examples  a coupled  rectan- 
gular and  circular  harbor  connected  to  the  open- sea.  Carrier,  Shaw 
and  Miyata  (1971)  also  theoretically  studied  resonant  oscillations  in 
harbors  of  special  shapes  (rectangular,  circular,  or  circular  sectors) 
connected  to  the  open-sea  through  a rectangular  entrance  channel  and 
found  that  for  inviscid  conditions  the  length  of  the  entrance  channel  affects 
both  the  wave  amplification  at  resonance  within  the  harbor  and  the 
frequency  bandwidth  of  the  response  near  resonance. 

Attempts  to  study  the  problem  of  wave  induced  oscillations  in 
harbors  of  complicated  geometry  began  with  the  two  independent  hydraulic 
model  studies  by  Knapp  and  Vanoni  (1945)  and  Wilson  (for  work  in  the 
period  of  1942  to  1951,  see  Wilson  (1959,  I960)).  Knapp  and  Vanoni's 
study  was  in  connection  with  harbor  improvements  at  the  Naval  Operating 
Base,  Terminal  Island,  California  (the  present  East  and  West  Basins  of 
Long  Beach  Harbor).  One  purpose  of  the  study  was  to  choose  the  optimum 
mole  alignment  and  in  this  connection  the  characteristic  response  of  the 

basin  was  also  determined.  Wilson's  work  dealt  with  the  problem  of 

« 

surging  in  Table  Bay  Harbor,  Capetown,  South  Africa  and  the  interaction 
of  moored  vesseslwith  wave-induced  oscillations  in  the  harbor.  Through 
experimental  studies  certain  proposed  modifications  were  suggested  for 
reducing  the  harbor  surging. 
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Although  hydraulic  model  studies  usually  can  provide  reliable 
information  on  harbor  surging  and  the  investigator  can  develop  corrective 
procedures  using  the  model,  generally  these  models  are  expensive  and 
require  a considerable  amount  of  time  to  operate  intelligently.  There- 
fore, many  researchers  have  been  searching  for  a method  to  theoretically 
determine  the  wave -induced  oscillations  in  harbors  of  arbitrary  shape. 

Such  theoretical  results  at  the  very  least  provide  a guide  for  initial 
calculations  and  for  the  experimental  program  if  indeed  a model  study 
is  deemed  necessary. 

Wilson,  Hendrickson,  and  Kilmer  (1965)  and  Leendertse  (1967) 
have  studied  numerically  the  three-dimensional  oscillations  in  bays  or 
harbors  of  arbitrary  shape  and  variable  depth  by  using  finite -difference 
techniques.  Both  methods  require  that  the  wave  condition  at  the  open 
boundary  (or  at  the  harbor  entrance)  be  assumed  or  specified  from  field 
data.  Recently  two  independent  studies  on  wave  oscillations  in  an  arbitrary 
shape  harbor  with  constant  depth  have  been  reported  by  Hwang  and  Tuck 
(1970)  and  Lee  (1969,  1970)  and  also  Lee  and  Raichlen  (1970).  Hwang 
and  Tuck  obtained  their  solution  by  superimposing  the  scattered  wave 
pattern  along  the  entire  reflecting  boundary  (including  the  coastline)  to 
the  standing  wave  system;  they  have  confirmed  their  analysis  by  com- 
paring with  the  experimental  results  obtained  by  Ippen  and  Goda  (1963) 
for  a fully  open  rectangular  harbor.  In  Lee  (1969)  the  domain  of  interest 
was  divided  into  two  regions:  a region  which  defines  the  limit  of  the 

harbor  and  the  open- sea  region.  Solutions  in  each  region  were  obtained 
in  terms  of  the  unknown  boundary  condition  at  the  entrance  with  the 
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response  determined  by  matching  the  solutions  at  the  harbor  entrance. 
This  theoretical  analysis  has  been  confirmed  experimentally  applying 
the  theory  to  two  circular  harbors,  a rectangular  harbor,  and  a constant 
depth  model  of  the  East  and  West  Basins  of  the  Long  Beach  Harbor. 


1.  2 DESCRIPTION  OF  PRESENT  STUDY 

The  present  study  is  an  extension  of  the  studies  by  Lee  (1969)  and 
Lee  and  Raichlen  (1970).  As  just  mentioned,  in  those  studies,  the  entire 
domain  was  divided  (at  the  harbor  entrance)  into  two  regions:  (1)  the 
harbor  region  which  forms  the  interior  limit  of  the  harbor  boundary  and 
(2)  the  open-sea  region.  The  solution  in  each  region  was  formulated  as 
an  integral  equation  in  terms  of  the  unknown  boundary  condition  at  the 
harbor  entrance.  A method  was  used  to  solve  for  the  integral  equation 
by  converting  it  to  a matrix  equation.  In  this  method  the  harbor  boundary 
was  divided  into  a sufficiently  large  number  of  straightline  segments, 
thus,  continuous  integration  in  the  solution  is  replaced  by  discrete 
summation  along  the  harbor  boundary.  The  final  solution  was  obtained 
by  matching  the  solutions  obtained  from  both  regions  at  the  harbor 
entrance.  This  method  has  been  successfully  applied  to  harbors  of 
several  geometrical  shapes  including  a complicated  harbor  geometry. 

In  principle,  this  method  can  be  applied  to  a harbor  of  any  arbitrary 
shape  connected  to  the  open- sea  as  long  as  the  water  depth  in  both 
regions  can  be  assumed  constant.  However,  in  order  to  realize  a 
prescribed  accuracy  for  the  numerical  calculations,  the  ratio  of  the 
length  of  the  straightline  boundary  segments  to  the  wave  length  must 
be  kept  within  certain  limits.  It  was  found  that  to  insure  good  experi- 
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mental  agreement  this  ratio  should  be  less  than  one-tenth.  This  require- 
ment implies  that  for  a harbor  with  a very  complicated  shape  or  short 
waves  the  storage  space  of  a digital  computer  must  be  quite  large  to 
handle  the  problem  and  the  computation  time  on  such  a computer  may  be 
fairly  long. 

In  order  to  improve  this  aspect  of  the  analysis  the  present  study 
was  initiated.  A method  has  been  developed  where  the  harbor  region  is 
divided  into  several  subregions,  and  the  same  method  as  just  described 
is  used  to  form  a solution  in  each  region  in  terms  of  the  unknown  boundary 
condition  at  the  common  boundary  which  separates  adjacent  regions.  In 
this  way  both  the  required  computer  storage  and  the  computation  time 
are  reduced  significantly.  As  an  example,  the  size  of  the  computer 
storage  for  the  computer  program  presented  by  Lee  (1969)  for  the  Long 
Beach  Harbor  Model  was  206,000  bytes  on  the  IBM  360/75.  (Each  eight 
binary  "bit"  is  referred  to  as  a byte;  four  bytes  form  a word.)  The 
execution  time  for  one  wave  number  (k  = 2.  35  ft  *)  including  the  com- 
pilation of  the  main  program  was  28.  9 sec.  Excluding  the  time  required 
for  compiling  the  main  program,  the  computation  time  was  24.  0 sec. 

By  applying  the  coupled-basins  method  to  compute  the  response  character- 
istics of  the  same  Long  Beach  Harbor  Model  (with  the  interior  harbor 
now  divided  into  two  regions:  the  East  Basin  and  the  West  Basin)  the 
computer  storage  on  the  same  IBM  360/75  computer  was  reduced  to 
145,000  bytes  while  the  execution  time  for  the  same  wave  number  was 
now  23.  4 sec  (including  the  compilation  of  the  main  program).  Sub- 
tracting the  time  for  compilation  of  the  program  from  this,  the  actual 
computation  time  was  13.0  sec.  Thus,  immediate  advantages  of  the 
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present  method  are  clearly  seen  and  more  complicated  harbor  geometries 
can  be  studied  for  a given  computer  size  as  well  as  investigating  shorter 
waves  than  could  be  investigated  with  the  method  of  Lee  (1969). 

The  coupled-basins  theory  which  was  developed  is  presented  in 
Chapter  2.  The  experimental  results  obtained  for  a model  of  the  East  and 
West  Basins  of  Long  Beach  Harbor  reported  by  Lee  (1969)  were  used  as 
initial  confirmation  of  the  present  theory.  In  addition,  experiments  were 
conducted  using  an  arrangement  of  two  connected  circular  basins  and  a 
harbor  which  consists  of  a circular  segment  with  a rectangular  channel 
connecting  it  to  the  open-sea.  In  Chapter  4 theoretical  results  are 
compared  to  these  experiments  for  the  three  cases. 

The  second  objective  of  this  investigation  was  to  study  in  detail 
certain  aspects  of  the  response  of  coupled  basins  of  simple  geometry 
compared  to  the  response  of  the  individual  basins  alone.  This  aspect  of 
harbor  resonance  is  quite  important  considering  that  one  method  of 
expanding  harbors  is  to  add  additional  slips  and  side  basins.  Before 
construction  takes  place  it  is  important  that  the  effect  of  these  basins  on 
the  main  harbor  be  completely  understood.  Several  different  arrangements 
of  circular -coupled  basins  and  circular-rectangular-coupled  basins  have 
been  investigated  theoretically  in  this  connection  and  the  analytical  results 
are  presented  in  Sections  4.  2 and  4.  3.  In  addition  to  the  response 
characteristics  of  such  harbors,  the  amplitude  distribution  in  a coupled- 
circular  harbor  is  presented  for  a particular  arrangement  for  several 
values  of  the  wave  number  parameter  ka  (the  product  of  the  wave  number 
k and  the  radius  a).  The  effect  of  a rectangular  entrance  canal  on  the 
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response  of  a circular  harbor  is  presented  in  Section  4.  3.  2.  This 
aspect  of  the  problem  was  not  pursued  in  detail,  but  was  investigated 
primarily  to  determine  the  effect  on  the  response  of  one  basin  with  a 
particular  shape  of  one  with  a different  shape;  these  results  indicate  a 
fruitful  area  of  further  study.  In  connection  with  this  preliminary  study 
of  entrance  channel  effects,  the  effect  of  the  side-channel  resonators 
on  the  response  of  harbors  with  entrance  channels  was  briefly  investi- 
gated theoretically  with  results  also  presented  in  Section  4.  3.  2.  A 
qualitative  discussion  of  the  effect  of  energy  dissipation  on  the  harbor 
response  is  presented  in  Section  4.  4 with  some  attention  given  to  the 
problem  of  scaling  model  results  to  the  corresponding  prototype. 
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CHAPTER  2 


THEORETICAL  ANALYSIS 


The  theoretical  analysis  for  the  wave  induced  oscillations  in 
harbors  which  are  composed  of  coupled  basins  with  arbitrary  shape  and 
a constant  depth  is  presented  in  this  chapter.  The  complete  harbor  is 
divided  into  several  regions,  and  the  solution  to  the  boundary  value  pro- 
blem in  each  region  is  formulated  as  an  integral  equation.  This  integral 
equation  is  approximated  by  a matrix  equation  which  can  be  solved  using 
a high-speed  digital  computer.  The  final  solution  is  obtained  using 
matching  conditions  at  the  harbor  entrance  and  the  common  boundaries 
between  regions,  i.  e.  , equating,  at  the  entrance  as  well  as  at  the 
common  boundaries,  the  wave  amplitudes  and  their  normal  derivatives 
obtained  from  the  solutions  in  each  region. 


2.  1 DEVELOPMENT  OF  THE  HELMHOLTZ  EQUATION 

Assuming  an  irrotational  flow,  a velocity  potential,  $,  can  be 
defined  such  that  the  fluid  particle  velocity  can  be  expressed  as  u = V6. 
Thus,  from  the  continuity  equation  for  an  incompressible  fluid,  Laplace's 


equation  is  obtained: 


V.u  = v $ = 0. 


A solution  of  $ is  sought  in  the  following  form: 


$(x,  y,  z;t) 


-^f(x,y)Z(z)e'AXJt, 


where  a is  the  angular  frequency  defined  as  2ir/T  (T  is  the  wave  period), 


*53? 


, , 


r--  _ | 
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X = y^T,  and  f(x,y),  termed  the  wave  function,  describes  the  variation 
of  $ in  the  x,y  direction.  (The  coordinate  system  is  defined  in  Fig.  1.  ) 


Substituting  Eq.  2 into  Eq.  1 one  obtains: 

r rv2  f \ i ry2  , 


i (b*i  ,9'f^  . i a"z 
f \3xa  + a7v  " zaP" 


(3) 


If  Eq.  3 is  equated  to  a constant,  say  -k3 , then  the  following  equations 
are  obtained: 


dsZ 


CT  - k‘Z  = 0 . 


dz‘ 


(4) 


9sf  , a2f  . .2,  n 
aja  + 5yS-  + k f = 0 • 


(5) 


The  boundary  condition  at  the  bottom  and  the  linearized  dynamic 
free  surface  condition  are  respectively: 


a$ 


ff(x,y,  -h:t)  = 0 , 


-jnt  1 f 

ri(x,y:t)  = A.f(x,y)e  ^ 


(6.  a) 
(6.b) 


z = 0 


where  the  depth  h is  assumed  constant,  q is  the  displacement  of  water 
surface  from  the  still  water  level,  is  the  amplitude  at  the  crest  of  the 
incident  wave,  and  g is  the  acceleration  of  gravity. 

The  function  Z(z)  which  satisfies  Eqs.  4 and  6 can  be  found  as: 


Z(z)  = 


A^g  cosh  k(h+z) 


cosh  kh 

Thus,  the  velocity  potential  $ becomes: 


$(x,y;t)  = jL  -1 


A.g  cosh  k(h+z) 


f(x,  y)  e 


-Xat 


1 


cosh  kh 


(7) 
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Substituting  Eqs.  6b  and  7 into  the  linearized  kinematic  free 
surface  condition,  8ri/8t  = (8®/8z)  _n,  the  well  known  "dispersion 

Z — U 

2 

relation"  for  water  waves  is  obtained:  0 = gk  tanh  (kh);  therefore, 

the  arbitrary  constant  k used  in  Eqs.  4 and  5 is  the  wave  number  defined 
as  2n7L  (L  is  the  wave  length). 

To  complete  the  expression  for  the  velocity  potential,  the  main 
problem  which  remains  is  to  determine  the  wave  function,  f(x,y),  which 
satisfies  Eq.  5,  known  as  the  Helmholtz  equation,  and  the  boundary 
condition  that  there  is  no  flow  through  solid  boundaries  (such  as  the 
coastline  and  the  boundary  of  the  harbor  in  this  problem)  and  also  the 
radiation  condition  which  will  be  discussed  later. 

2.  2 SOLUTION  OF  THE  HELMHOLTZ  EQUATION  FOR  COUPLED  BASINS 

The  procedure  for  determining  the  wave  function,  f(x,y),  and  thus 
the  response  of  the  harbor  to  periodic  incident  waves  can  be  outlined  as  follows: 

(i)  The  domain  of  interest  shown  in  Fig.  2 is  divided  into 

several  regions:  the  infinite  ocean  region  (Region  I), 

and  regions  bounded  by  the  limits  of  the  harbor  and  various 
interior  divisions  (Region  II- 1,  Region  II-2,  Region  II-3). 

(ii)  The  function  fT  in  Region  I is  expressed  in  terms  of  the 
unknown  value  of  the  normal  derivative  8fx  /8n  at  the 
harbor  entrance.  The  function  f21  in  Region  II- 1 is  expressed 
in  terms  of  the  unknown  value  of  8f2l  /8n  both  at  the  harbor 
entrance  and  at  the  common  boundaries  between  other  basins 
in  Region  II,  e.  g.  , line  CD  and  EF  in  Fig.  2.  Similarly,  the 
function  f22  in  Region  II-2  (or  the  function  f23  in  Region  U-3) 


Fig.  2 A Definition  Sketch  of  an  Arbitrary  Shape  Harbor  Containing 
Coupled  - Basins 
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is determined  in  terms  of  the  unknown  value  of  9f22  /9n  (or 
9f23/9n)  at  the  common  boundary  between  Regions  II- 1 
and  II-2  (or  between  Regions  II- 1 and  II-3).  It  is  noted 

that  the  functions  of  f*  , fsl  , f22 all  must  satisfy 

the  Helmholtz  equation,  ^ f + k f = 0,  in  the  respective 
region. 

(iii)  At  the  boundaries  between  regions  the  wave  functions  and 
the  normal  derivatives  are  equated  to  solve  ultimately  for 
the  value  of  the  normal  derivative  at  the  harbor  entrance, 

9fx  / 9n.  This  is  denoted  as  the  "continuity"  condition  and 
can  be  formulated  as: 

fx  = f 2 1 , 9fx  / 9n  = -9f21  /9n  at  the  entrance  to  the  open- 

sea  (boundary  AB  as  shown 
in  Fig.  2), 

f2 1 = f22,  9f21 /9n  = -9f22/9n  at  the  imaginary  boundary 

between  Regions  II- 1 and  II-2, 
(boundary  C D as  shown  in  Fig. 


2 ) , and 

f2i  = f23,  9f21  / 9n  = -9f23/9n  at  the  imaginary  boundary 

between  Regions  II- 1 and  11-3 
(boundary  EF  as  shown  in  Fig. 
2). 

This  matching  procedure  is  done  simultaneously. 

The  negative  sign  in  the  normal  derivative  appearing  in  the 
formulated  matching  conditions  is  the  result  of  the  adapted 
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(iv)  Once  the  normal  derivative,  8f/8n,at  each  common  boundary 
is  determined,  the  wave  function  f within  each  basin  can  be 
evaluated. 


This  method  can  be  applied  to  an  arbitrary  number  of  coupled 
interior  basins,  although  the  example  to  be  discussed  (Fig.  2)  only  has 
three  obvious  interior  regions. 


2.  2.  1 The  Function  f2i  in  Region  II- 1 

The  function  f21  at  any  position  x inside  Region  II- 1 can  be 
expressed  in  terms  of  the  values  of  f21  and  9f2i  /8n  at  the  boundary  of 
the  region  by  applying  Weber's  solution  of  the  Helmholtz  equation  (see 
Baker  and  Copson  (1950)  or  Lee  (1969)): 


s 

where  xq  is  the  position  vector  of  a point  on  the  boundary,  H^(kr)  is  the 
zero  order  Hankel  function  of  the  first  kind,  n is  directed  outward  and 
normal  to  the  boundary,  and  r is  the  distance  |x-xq|.  The  contour 
integration  is  to  be  performed  along  the  boundary  of  the  region  moving 
in  a counterclockwise  direction. 

“4 

Eq.  8 shows  that  the  wave  function  f21  (xq)  must  be  known  at  all 
boundaries  of  Region  II- 1 before  the  function  f21(x)  at  any  interior  point 
can  be  obtained.  To  ->btain  a solution,  the  wave  function  at  the  boundary, 
fa!  (xQ),  is  expressed  in  terms  of  the  normal  derivative  of  the  wave  function 
at  the  boundaries  between  Region  II- 1 and  Regions  H-2  and  II-3  and  that  at 
the  harbor  entrance.  This  is  accomplished  by  modifying  Eq.  8 by  allowing  the 
field  point  x to  approach  a boundary  point  xi(x.,y.)  from  the  interior  of 
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the  region  (see  Fig.  2).  Thus,  if  the  boundary  is  sectionally  smooth, 

Eq.  8 becomes  (see  Lee,  1969  for  this  derivation): 

tHo  l(kr,]  - Ho  • <9> 

s 

where  r = x.-x  j. 

1 o' 

This  integral  equation  (Eq.  9)  cannot  be  solved  analytically  for 
an  arbitrarily  shaped  boundary:  however,  an  approximate  method  can 
be  applied.  Eq.  9 is  expressed  in  discrete  form  by  dividing  the  boundary 
of  Region  II- 1 into  a sufficiently  large  number  of  straight  line  segments 
(Nx): 


Nt 


fSl  (Xj)  = 


P iAs. 

:4 


^ ^ g 

‘ 2h' 1 >(kr>dr. 

(10) 


where  r.j  = |x.-Xj|.  The  notation  and  the  numbering  system  used  for 
the  straight-line  segments  are  shown  in  Fig.  3.  The  segments  are 
numbered  starting  at  the  right-hand-corner  of  the  harbor  entrance  and 
then  proceeding  counterclockwise.  (It  should  be  noted  that  for  the  other 
two  regions.  Region  II-2  and  Region  II-3,  the  numbering  system  starts 
at  the  right-hand -corner  of  the  common  boundary  between  regions  and 
then  also  progresses  counterclockwise.  ) 

Eq.  10  can  be  expressed  in  matrix  form  as: 


Fx  = -f(GnFi  -GPx)  , 


(ID 


Note:  Entrance  to  the  open-sea:  p segments 

Common  boundary  between  Regions  II- 1 and  II-2:  di  segments 
Common  boundary  between  Regions  II- 1 and  II-3:  ds  segments 


Fig.  3 A Definition  Sketch  of  the  Harbor  Boundary  Approximated  by 
Straight-Line  Segments 


where: 


= fsifXj) 

(i  = 1.2. 

Ni), 

(12.  a) 

(i  = 1,2. 

Ni), 

(12. b) 

<Gn>ij  = 


y\2i 

r. . 

J iJ  J 

(for  i,  j = 1,2 


J 11  i ii  i** 

Nj  ; i^j),  (12.  c) 

(Gn'ii  * ;(U0  -0&K  «°ri  = >>2 N>1-  <12'd> 

(fori.j  = 1,2,  Mi  ; ii^j ),  (12.  e) 

(G)..  = fl  +x|[log(^i)-  0.42278]}as. 

(for  i = 1,2 Nj  ) . (12.  f) 


,G»ij  - «" ’<kv4*j 


,kAs. . 


For  the  derivation  of  the  diagonal  elements  of  the  matrices 
and  G (Eqs.  12.  d and  12.  f),  the  asymptotic  expressions  for  the  Hankel 
functions  are  used,  i.  e.  , H|^(kr)  ~ an<*  H^^kr)  ~ 1 

(log^rj-  + y)  as  kr  -*0  (y  = 0.577216,  Euler's  constant).  Detailed 
derivations  of  Eqs.  12  are  given  by  Lee  (1969). 

The  vector  Px  in  Eq.  11  involves  the  unknown  value  of  9f21  /9n 
at  the  harbor  entrance  and  at  the  common  boundaries  between  Regions 
II- 1 and  H-2,  and  Regions  II- 1 and  U-3.  Since  the  value  of  9f2i  /9n  at 
solid  boundary  is  zero,  the  vector  Px  can  be  expressed  as  follows: 
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Eq.  13  can  be  written  in  a simpler  fashion  as: 

Pi  = Ui  • C , (14) 

where  U*  is  a Nj  x D matrix  as  defined  in  Eq.  13,  (in  which  D = p + di  + d2 
and  p is  the  total  number  of  segments  at  the  entrance  to  the  open-sea,  dx  is 
the  total  number  of  segments  at  the  common  boundary  between  Regions 
II- 1 and  II-2,  and  d2  is  the  total  number  of  segments  at  the  common 
boundary  between  Regions  II- 1 and  II- 3).  The  vector  <3  in  Eq.  14  consists 
of  the  unknown  values  of  3f2l  /9n  at  the  harbor  entrance,  and  at  the  common 
boundaries  between  Regions  II- 1 and  II-2  as  well  as  between  Regions  II- 1 and 
U-3. 

Substituting  Eq.  14  into  Eq.  11  and  rearranging,  one  obtains: 

B - (iO„  + 1)’l(?ou,)c 

= MtC  , (15) 

where  M!  = +1  ) (^rGUi  ) is  a N]  x D matrix  and  can  be  computed 

directly,  since  the  matrices  G^,  G,  Ua  are  known  matrices  at  this  stage 
of  the  development  (I  is  an  identity  matrix). 

Eq.  15  shows  that  the  value  of  f21  (x^)  on  the  boundary  of  Region 
II-  1 can  be  expressed  as  a function  of  the  values  of  0f2  j / 3n  at  the  harbor 
entrance  and  those  at  the  common  boundaries.  Th»  vector  , i.  e.  , the 
normal  derivatives,  will  be  determined  through  a matching  procedure 
which  will  be  described  in  Subsection  2.  2.  5. 
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2.  2.  2 The  Function  f2a  in  Region  II -2 

Following  a procedure  which  is  similar  to  that  used  in  Sub- 
section 2.2.  1,  the  function  f22(x)  at  any  position  x inside  Region  II- 2 can 
be  expressed  in  terms  of  the  values  of  f22  and  8f22/3n  on  the  boundary 


4.  A = - Hol,kr|^£”<;;o|]ds(*o)  • 1161 

s 

where  r = |x  - x . 

' o 

By  allowing  the  interior  point  x to  approach  a boundary  point  x^ 
the  following  integral  equation  (analogous  to  Eq.  9)  is  obtained: 

WJi>  = -?J[£«<*0)^<1)(kr)-H»,(kr)^£ai,(:0)]d.(;0)  . (17) 

s 

where  r = |x.  - x I . 

1 l o' 

The  approximate  method  used  in  Region  II- 1 is  applied;  thus  Eq. 

17  can  be  written  as  a matrix  equation: 

F2  = -j  (GnFt  - GPp ) . (18) 


where: 


F2  = f2  2 (X.  ) 


P2  = 


i = 1,  2. 


i = 1,  2, 


N2  ; 


the  matrices  G and  G are  each  an  N2  x N2  matrix.  The  elements  of 
n 

the  matrices  G and  G can  be  calculated  from  the  same  expressions 
n 

given  in  Eqs.  12c,  12d,  12e,  12f.  Of  course,  it  is  realized  that  the  index 
i or  j referred  to  in  these  equations  represents  the  boundary  points  of 
the  Region  II-2. 


The  vector  Ps  involves  the  unknown  values  of  9f3S/8n  at  the 
common  boundary  between  Region  II- 1 and  Region  II-2.  Thus,  using 
the  conditions  described  at  the  beginning  of  Section  2.  2 (see  iii)  the 
vector  Ps  can  be  related  to  the  value  of  8f21  /8n  at  the  common 
boundary. 

The  vector  P2  can  thus  be  written  as  follows: 


where  the  matrix  U2  is  a N2  x dt  matrix  with  the  antidiagonal  elements 
of  the  first  dx  rows  and  dx  columns  equal  to  -1  and  with  the  other 


i = 1,  2.  .. 


where  F3 

P3 


fs^Xj)  , 
glWXj), 


i = 1.  2, 


n3 

N3 , and 


the  matrices  and  G are  both  N3  x N3  matrices,  their  elements  can  be 
calculated  by  applying  the  formulas  which  are  given  in  Eq.  12c,  12d, 

12e,  and  12f  to  the  boundary  points  of  Regions  II-3. 


where  the  matrix  U3  is  a N3  x ds  matrix  with  the  antidiagonal  elements 
of  the  first  ds  rows  and  d2  columns  equal  to  -1  and  the  remaining  elements 
equal  to  zero.  The  vector  C3  is  a ds  dimensional  vector  with  its  elements 
equal  to  the  last  da  elements  of  the  vector  C defined  in  Eq.  13. 


Substituting  Eq.  23  into  Eq.  22  and  solving  for  the  vector  F3,  one 


obtains: 


F3  = (£cn  + l)  (jGU3)£  = M3C3  . 


where  M3  is  a N3  x d2  matrix  which  can  be  calculated  directly. 

If  the  region  forming  the  boundary  of  the  harbor  is  divided  into 
more  than  three  regions,  i.  e.  , more  than  Region  II- 1 and  II-2,  and  II-3, 
an  approach  similar  to  that  used  in  Subsections  2.  2.  2 and  2.  2.  3 can  be 
used  to  formulate  the  solutions  for  additional  regions. 

In  order  to  solve  for  the  unknown  vector  C defined  in  Eq.  13,  the 
solution  in  Region  I (the  open-sea  region)  must  first  be  obtained.  This 
is  presented  in  the  following  subsection. 

2.2.4  The  Function  fx  in  Region  I (Open-Sea) 

The  solution  fj  of  the  Helmholtz  equation  in  Region  I can 
be  considered  as  composed  of  three  separate  parts:  a function  f. 

representing  an  incident  wave,  a function  f representing  a reflected 
wave  considered  to  occur  as  if  the  harbor  entrance  were  closed,  and  a 
function  f3  (termed  the  radiated  wave  function)  representing  a correction 
to  f due  to  the  presence  of  the  harbor. 

Thus  the  wave  amplitude,  r^fx.yjt)  in  Region  I can  be  expressed 


r)  i (x,  y;t)  = A.  fj  e 


-a.  at 


= A.(f.  + f + f )e 
li  r s' 


-x,at 


The  incident  wave  function,  f^,  can  be  specified  in  an  arbitrary 
fashion;  for  example,  a periodic  incident  wave  with  the  wave  ray  at  an 
angle  a to  the  coastline  (x-axis  as  shown  in  Fig.  2)  can  be  represented 
as  f.  (x,y)  = exp[,tk(x  cos  a + y sin  a)].  The  reflected  wave  function  f 


can  be  obtained  from 


incident  wave  with  wave  ray  perpendicular  to  the  coastline  (a  = 90  ),  the 
function  f. (x,y)  can  be  represented  by  £e^^  (the  factor  £ is  chosen  for 


convenience).  In  the  following  discussion,  the  incident  waves  will  be 


considered  as  normally  incident  to  the  coast  line;  this  condition  also  was 


treated  experimentally  in  this  study. 

The  major  problem  in  defining  the  open- sea  wave  system,  i.  e 


function  f,  , is  to  evaluate  the  radiated  wave  function  f,  which  must  also 


satisfy  the  Helmholtz  equation 


and  the  following  boundary  conditions 


at  the  coastline  (boundary  AO  and  BO'  of  Fig.  2) 


si  at  the  harbor  entrance  (boundary  AB  of  Fig.  2),  ai 


this  is  called  the  radiation  condition 


(where  r 


The  function  f3  which  satisfies  Eq.  26  and  the  boundary  condition  is 
obtained  using  Weber's  solution  of  the  Helmholtz  equation: 
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where  xq  is  the  source  point  (xq,0)  along  the  x-axis,  x is  the  field  point 
(x,y)  in  the  open-sea  region,  r = ~ xQ)a  + ya  , and  the  integration  is 

to  be  performed  along  the  x-axis.  It  should  be  noted  that  the  fundamental 
solution  H^(kr)  is  necessary  to  satisfy  boundary  condition  (iii)  above, 
i.  e.  , the  radiation  condition,  (see  Lee,  1969  for  a detailed  discussion). 

In  order  to  solve  for  the  value  of  f3  at  the  harbor  entrance  in 
terms  of  the  value  of  9f3/9n  at  the  entrance,  the  field  point  x is  allowed  to 
approach  the  x-axis  at  the  point  (x.,  0),  thus,  from  Eq.  27  one  obtains  the 
following  integral  equation  (similar  to  Eq.  9): 

f3(x.,0)  = J[f3Uo,0)^Hj)1)(kr)  - H^)(kr)^f3(xo.0)]ds(xo,0)  , (28) 

s 

where  r = lx.  - x . 

1 i o' 

Substituting  the  boundary  conditions  associated  with  Eq.  26  into 
Eq.  28,  and  expanding  the  terms  inside  the  integral,  the  following  simplified 
equation  results: 


f3(xi,0)  = 


y f H(1)(k|x.  -x  I)#  [f21(x  ,0)]dx 
2 J o ' l o 1 9n  L 2 1 o o 


AB 


« -T  7 H<U(k|x.  - x | ) C.  As . -yC.  r^^si2ri  + ^(log(^)  + y)]dr  , 

£ O 1 1 1 1 L 1 o TT  C, 

. J J J Q 

(29) 

th  tli 

where  x^  and  x^  are  the  mid -points  of  the  i and  j segments  of  the  harbor 
entrance  respectively,  As^  is  the  length  of  the  j1*1  segment,  the  term  C.  in 
Eq.  29  is  the  value  of  9f2l  /9n  at  the  mid-point  of  the  j**1  entrance  segment, 


and  p is,  as  mentioned  before,  the  total  number  of  segments  into  which 
the  entrance  is  divided.  (The  symbol  y in  Eq.  29  is  referred  to  as 
Euler's  constant  which  is  equal  to  0.  577216.  ) 

Once  the  value  of  f3(x.,0)  at  the  entrance  is  determined  from 
Eq.  29,  the  value  of  f1(x^,0)  at  the  harbor  entrance  can  be  obtained  by 
adding  the  contribution  of  f.  and  f to  f3 ; thus: 

P 

fa(x.)  = 1 + 7 H..  C.  (30) 

j=l 

for  i = 1,2,  ...  p.  The  first  term  on  the  right-hand  side  of  Eq.  30 

represents  the  incident  plus  reflected  wave  if  the  entrance  is  closed. 

This  amplitude  is  unity  since  for  convenience  the  amplitude  of  the 

incident  wave  was  chosen  as  one-half.  The  matrix  H..  in  Eq.  30 

represents  the  contribution  of  the  function  f3  and  it  is  rewritten  in 

matrix  form  from  Eq.  29,  where  H..  = - \kr. . ) As ; for  i,j  = 1,2, 

n lj  2 o ij  J J 

...p;i  4 j,  and  H..  = - j[  1 + ^log'^^p )-  0.  42278  )~*Aa,  for  i = 1,2, 

• • • P* 

At  this  stage  of  the  development  the  solution  of  the  Helmholtz 
equation:  V f + k f = 0,  for  each  region  has  been  formulated  in  terms 
of  the  value  of  the  unknown  derivative  af/9n  at  each  common  boundary 
between  interior  regions  and  at  the  harbor  entrance.  In  the  next  subsection, 
a matching  procedure  is  discussed  which  results  in  a solution  for  the 
harbor  response  for  the  individual  basins  from  these  solutions. 

2.  2.  5 Matching  Solutions  at  the  Harbor  Entrance  and  the  Common 
Boundaries  Between  Regions 

The  solutions  for  the  wave  function,  f,  for  the  various 


regions  are  summarized  as  follows: 


in  Eq.  3 If  are  interchanged  so  that  they  will  correspond  to  those  of 
Eq.  31e. 

The  matching  procedure  consists  of  equating  Eq.  31a  to  Eq.  31b, 
Eq.  31c  to  Eq.  31d,  and  Eq.  31e  to  Eq.  3 1 f simultaneously.  Thus,  one 
obtains  the  following  matrix  equation: 


MAc  = s + HAc 


(32) 


where: 


1)  Ma  is  a D x D matrix  consisting  the  first  p rowt,  the  nij 


th 


row  to  na*k  row,  and  the  m2t^1  row  to  na^  row  of  the  matrix  whose 


/ 


ma  - 


are  defined  in  Eqs. 

12,  14,  and  15: 

(Mi  j (Ma  )12 

• • • • 

• • • • • • 

• • • • 
(Mi  )pl  (Mi  )p2  . . 

• • • • • • 
<m>>pd 

(M^mil(Ml)mi2-  ’ 

'^'miD 

(M^nil  <M*>nl2  * * 

(M^nlD 

<Mi>m2l<M>>m22-  ' 
• • • • 



(Mj)  (Mi)  . . 

rig  1 ns  u 

(Ml  ^ D 

\ 


2)  C is  a D x 1 column  vector  consisting  ot  the  values  of  3f21  /9n  at 


every  segment  representing  the  harbor  entrance  and  the  common  boundaries 
between  regions  and  is  defined  as: 


- 33  - 


I Cg  . . c c • c c 

L p mi  r»i  mg 


. C 

n®  J 


3)  S is  a D x 1 column  vector  with  its  first  p elements  repre- 
senting the  value  of  the  incident  wave  plus  the  reflected  wave  at  the 
harbor  entrance,  these  are  equal  to  unity,  and  with  the  other  elements 
equal  to  zero.  Thus,  it  is  defined  as:  [11..10...0] 

4)  is  a D x D matrix  defined  as: 


Hi,  H, 


HlZ 

Hip  0 

0 o 

0 

H22 

h2p  . 

. 

• • 

• 

• • 

. 

• • 

H „ 

H 0 

0 0 

0 

P2 

PP 

°(M,  >dil 

• <M=>dldl 0 

0 

• • 

• • • 

• 

0 (Ma)n 

• <“.)«,  0 

0 

0 0 

0 <M3>dal 

• ,M’»dad 

• • 

• • • 

• • 

0 0 

0 (M3)n 

• IM=!lds 

Thus,  Eq.  32  can  be  solved  for  the  unknown  vector  as  follows: 

, v-1 


£ * (ma  ‘ ha)  • 2 


The  unknown  value  of  9fal  /9n  at  the  harbor  entrance  is  equal  to 
the  first  p elements  of  the  vector  C,  the  value  of  9fal  /9n  at  the  common 
boundary  between  Regions  II- 1 and  II- 2 is  equal  to  the  next  dj  elements 
of  the  vector  C,  and  the  value  of  9fsl  /9n  at  the  common  boundary  between 
Regions  II- 1 and  U-3  is  equal  to  the  last  ds  elements  of  the  vector 


.. 
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Once  the  vector  £ is  obtained  from  Eq.  33,  the  value  of  the 
wave  function  at  the  boundary  of  each  region  can  be  determined.  The 
wave  function  f21  at  the  boundary  of  Region  II- 1 can  be  obtained  from 
Eq.  15,  while  the  value  of  f22  at  the  boundary  of  Region  II-2  can  be 
obtained  from  Eq.  20  and  f23  at  the  boundary  of  Region  II- 3 can  be 
obtained  from  Eq.  24. 

The  value  of  wave  function  f at  any  position  within  each  region 

can  be  determined  after  the  vector  C and  the  wave  function  f at  the 

boundary  of  each  region  are  evaluated.  For  example,  for  an  arbitrarily 

chosen  point  x in  Region  II- 1 the  value  of  f21(x)  can  be  calculated  by  the 

following  discrete  form  of  Eq.  8: 

Nj  D 

UiM  = -ff£  ^(x.jf-kH^Wjf^As.  -^H^WlC-As.}  (34) 

j=l  “ j=l 

where  x^  is  the  mid -point  of  the  j^1  boundary  segment  in  Region  II- 1,  and 
r = |x  - Xj  |.  In  an  analogous  manner,  similar  expressions  can  be  used 
to  determine  the  wave  functions  in  Regions  II-2  and  II- 3. 

The  response  of  a harbor  to  incident  waves  is  described  by  a 
parameter  called  the  "amplification  factor".  This  is  the  ratio  of  the 
wave  amplitude  T]2i(x,y)  at  any  position  (x,y)  inside  the  harbor  (i  = 1 
refers  to  Region  II- 1,  i = 2 refers  to  Region  II-2,  and  i = 3 refers  to 
Region  II-3),  to  the  sum  of  the  incident  plus  the  reflected  wave  amplitude 
at  the  coastline  with  the  harbor  entrance  closed: 


|Aifsi(x,y)e-'Wt| 

R yxsr,  * . . .^t,  = l£»i (*-y>l  (35) 


lAjIfj  + fr>.  '|  |Ai-I-.-tot| 


CHAPTER  3 


EXPERIMENTAL  APPARATUS  & PROCEDURES 

In  general  the  experimental  equipment  used  for  this  study  is 
described  by  Lee  (1969)  and  only  will  be  summarized  here.  An  over- 
all view  of  the  experimental  equipment  is  shown  in  Fig.  4. 

The  wave  basin  used  is  1 ft  9 in.  deep,  15  ft  5 in.  wide,  and 
31  ft  5 in.  long  with  vertical  walls  of  3/4  in.  marine  plywood  and  a 
floor  constructed  of  1 in.  marine  plywood.  The  bottom  of  the  wave 
basin  is  horizontal  to  within  at  least  ±0.  02  in.  and  has  been  treated  with 
a layer  of  polyester  resin  approximately  1/8  in.  thick  to  provide  a 
water-tight  seal. 

The  wave  generator  is  of  the  pendulum  type  1 1 ft  8 in.  long  and 
2 ft  high  located  at  one  end  of  the  basin,  and  is  designed  to  operate  either 
as  a paddle-  or  piston-type  wave  generator;  for  a detailed  description, 
the  interested  reader  is  referred  to  Raichlen  (1965).  The  wave  gener- 
ator is  driven  by  two  arms  connected  to  independent  cranks  which  in 
turn  are  connected  through  a pulley  system  to  a 1-1/2  hp  variable  speed 
motor.  The  cranks  allow  for  a maximum  stroke  of  ±6  in.  and  they 
can  be  adjusted  to  within  0.  001  in.  of  each  other.  Wave  periods  ranging 
from  0.  34  sec  to  3.  8 sec  can  be  obtained  with  this  system. 

The  wave  period  is  determined  by  measuring  the  rotational  speed 
of  one  crank  of  the  wave  generator.  This  is  accomplished  by  attaching 
a disc  with  360  evenly  spaced  holes  arranged  at  its  outer  edge  to  the 


Fig.  4 An  Over-All  View  of  the  Experimental  Equipment 
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water  maintained  between  6.  2 to  6.  5 to  insure  its  proper  functioning  as 
a corrosion  inhibitor. 

The  wave  filter  and  wave  absorbers  used  have  reflection  coefficients 
which  are  less  than  20%  for  the  experimental  range  of  depths,  wave  heights, 
and  wave  periods.  This  means  that  the  amplitude  of  a wave  which  passes 
through  the  absorber  (or  filter)  and  reflects  from  the  wall  (or  wave 
machine)  and  then  passes  back  through  the  absorber  (or  filter)  to  the 
main  wave  basin  is  reduced  by  more  than  80%.  A more  detailed  discussion 
of  the  characteristics  of  the  dissipation  system  is  given  by  Lee  (1969). 

Experiments  were  conducted  in  a well  defined  coupled  harbor  to 
test  the  theory  which  is  described  in  Chapter  2.  The  harbor  consisted  of 
two  circular  basins  with  the  same  diameter  connected  together  on  the  axis 
which  is  perpendicular  to  the  "coastline”.  The  first  basin  of  the  harbor 
had  a 60°  opening  which  communicated  directly  with  the  "open- sea"  and 
a 10°  opening  diametrically  opposite  connected  to  the  10°  opening  of  the 
second  basin  of  the  harbor.  Each  harbor  was  1 ft  6 in.  in  diameter  and 
1 ft  3 in.  high.  The  model  used  is  shown  in  the  photograph  presented  in 
Fig.  5.  It  is  seen  in  Fig.  5 that  the  two  cylinders  which  make  up  the 
harbor  are  connected  top  and  bottom  to  lucite  reinforcing  plates  which 
were  used  to  keep  the  basins  circular.  The  edges  of  the  entrance  between 
the  two  basins  and  at  the  "open-sea"  do  not  have  a zero  thickness  as 
assumed  in  the  theory,  but  due  to  construction  limitations  this  thickness 
was  approximately  1/2  in.  for  the  former  and  1/4  in.  for  the  latter.  The 
two  vertical  lucite  plates  to  either  side  of  the  main  harbor  entrance  shown 
in  Fig.  5 fitted  into  the  existing  end  wall  of  the  basin  to  form  a perfectly 


L 


f 


I 


<!> 


CHAPTER  4 

PRESENTATION  AND  DISCUSSION  OF  RESULTS 

Theoretical  results  obtained  using  the  method  developed  in 
Chapter  2 are  presented  here  along  with  some  limited  experimental 
results.  The  effect  of  interconnected  basins  in  a harbor  on  the  overall 
harbor  response  is  discussed  relative  to  the  response  of  the  individual 
basins. 

4.  1 A MODEL  OF  THE  EAST  AND  WEST  BASINS  OF  LONG  BEACH 
HARBOR 

A model  of  the  East  and  West  Basins  of  Long  Beach  Harbor  was 
studied  previously,  both  theoretically  and  experimentally,  by  Lee  (1969). 
The  theoretical  results  of  that  study  were  obtained  by  dividing  the  domain 
of  interest  into  two  regions:  the  open-sea  region  (Region  I)  and  the 
harbor  domain  (Region  II).  The  boundary  of  Region  II  was  divided  into 
75  segments  (with  two  segments  for  the  harbor  entrance)  to  obtain  a 
numerical  solution;  this  subdivision  is  shown  in  Fig.  6.  (The  theory 
presented  in  Lee  (1969)  was  referred  to  therein  as  the  "arbitrary-shape 
harbor  theory";  herein  it  will  be  referred  to  as  the  "single-basin 
theory".  ) 

The  validity  and  accuracy  of  the  coupled -basins  theory,  presented 


in  Chapter  2,  was  tested  initially  using  the  model  of  Long  Beach  Harbor 
which,  as  can  be  seen  in  Fig.  6,  is  essentially  composed  of  two  basins 
(the  East  and  the  West  Basins).  For  the  analytical  model  the  harbor  was 


N=75  (TOTAL  NO.  OF  SEGMENTS) 


Configuration  Used  in  Single  Basin  Theory  for  Analytical  Model  of 
East  and  West  Basins  of  Long  Beach  Harbor  California 


divided  at  the  junction  of  the  East  and  West  Basins;  this  division  is 
shown  in  Fig.  7.  The  boundary  of  Region  II- 1 (the  West  Basin)  is 
divided  into  50  segments  with  two  segments  for  the  entrance  to  the  open- 
sea,  and  four  segments  for  the  common  boundary  between  the  East  and 
West  Basins.  The  boundary  of  Region  II-2  (the  East  Basin),  shown  in 
Fig.  7,  is  divided  into  34  segments  including  four  segments  for  the 

i 

common  boundary  between  Regions  II- 1 and  H-2. 

Response  curves  have  been  obtained  using  the  coupled-basins 
theory  and  are  presented  in  Figs.  8 through  11  for  four  different  locations 
in  Long  Beach  Harbor.  The  ordinate  is  the  ratio  of  the  wave  amplitude 
at  the  particular  location  to  the  open- sea  standing  wave  amplitude  (Eq. 

35)  and  the  abscissa  is  the  product  of  the  wave  number  and  a characteristic 
harbor  length,  a,  (in  the  model  a = 1.44  ft,  in  the  prototype  a = 6768  ft). 

In  these  figures,  both  the  experimental  results  and  the  theoretical  results 
for  the  harbor  treated  as  one  basin  are  presented  for  comparison  (see 
Lee  (1969)).  Figs.  8 through  11  show  that  the  present  results  agree  well 
with  the  previous  single-basin  theory  throughout  the  full  range  of  ka  which 
was  investigated.  Only  small  differences  between  the  theories  can  be 
seen  in  the  response  at  resonance  which  is  the  region  where  both  theories 
are  somewhat  unreliable  because  both  theories  neglect  viscous  and  non- 
linear effects.  Actually  the  maximum  difference  between  the  two  theoretical 
results  is  less  than  3%.  Since  both  theoretical  curves  agree  well  with  each 
other  and  with  the  experimental  results  at  all  four  locations,  it  can  be  con- 
cluded that  the  wave  amplitude  distribution  within  the  harbor  is  also  pre- 
dicted correctly  by  the  coupled -basins  theory. 


0.2  ft. 


Fig.  7 Configuration  Used  in  Coupled- Basins  Theory  for  Analytical  Model  of 
East  and  West  Basins  of  Long  Beach  Harbor,  California 


Fig.  8 Response  Curve  at  Position  A of  the  Long  Beach  Harbor  Model 


Fig.  9 Response  Curve  at  Position  B of  the  Long  Beach  Harbor  Model 


Fig.  10  Response  Curve  at  Position  C of  the  Long  Beach  Harbor  Model 
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Fig.  1 1 Response  Curve  at  Position  D of  the  Long  Beach  Harbor  Model 
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4.  2,  THE  RESPONSE  OF  A CIRCULAR  BASIN  WITH  A 10°  OPENING 
CONNECTED  TO  A CIRCULAR  HARBOR  WITH  A 60°  OPENING 

Theoretical  and  experimental  results  have  been  obtained  for  a 
coupled -harbor  consisting  of  two  circular  basins:  the  first  with  a 10° 
opening  connected  to  a second  basin  which  has  a 60°  entrance  which 
opens  directly  on  the  open-sea  (see  Fig.  5).  The  response  curves  for 
two  positions  in  each  of  the  basins  are  presented  in  Figs.  12  through  15 
where  the  ordinate  is  the  amplification  factor  and  the  abscissa  is  the 
wave  number  parameter  ka  (the  product  of  the  wave  number  and  the 
harbor  radius).  In  the  experiments  the  radius  of  each  circular  basin 
is  0.  75  ft  and  the  depth  of  the  water  in  both  the  harbor  and  the  "open- 
sea"  was  1.0  ft. 

The  boundary  of  Basin  A was  subdivided,  in  applying  the  coupled- 
basins  theory,  into  38  segments  with  six  of  these  at  the  harbor  entrance 
and  two  at  the  common  boundary  between  Basins  A and  B.  Basin  B was 
divided  into  37  segments  including  two  segments  at  the  boundary  between 
Basins  A and  B. 

Figs.  12  through  15  show  that  the  theoretical  results  agree  reason- 
ably well  with  the  experiments  except  near  resonance  for  the  range  of  ka 
covered.  Near  the  various  resonant  modes  of  oscillation  the  amplification 
predicted  by  theory  is  greater  than  that  measured  experimentally;  how- 
ever, the  resonant  frequencies  are  accurately  predicted  by  the  theory. 

The  differences  in  amplification  mentioned  are  probably  due  to  viscous 
dissipation  which  affects  the  experiments  but  does  not  enter  the  inviscid 


A 


Fig.  12  Response  Curve  at  Center  of  Basin  A for  Two  Coupled-Circular  Harbors 


v 


coupled-basins  theory.  Similar  differences  were  reported  by  Lee  (1969)  in 


connection  with  '.he  application  of  the  single -basin  theory.  It  is  interesting 
to  note  that  some  of  the  largest  responses  are  associated  with  peaks  with 
very  narrow  bandwidths,  and  in  the  experiments  these  responses  are 


significantly  attenuated.  In  fact,  in  some  cases  these  maxima  have  almost 


Thus 


can  play  an  important  role  in  minimizing  resonance  effects.  Some 


comments  on  the  effect  of  viscous  dissipation  on  harbor  response  will  be 


ho  we  ver 


some  of  these  resonant  conditions  are  not  immediately  apparent  from  the 


figures.  The  reason  for  this  is  that  Figs.  12  through  15  refer  to  the 


response  at  a particular  location  in  one  harbor  as  a function  of  the  wave 


number  parameter  ka.  Hence  the  effect  of  resonance  for  a particular 


wave  number  at  one  position  in  one  basin  may  be  important,  but  at 


another  location  no  effect  may  be  seen.  Consider  Fig.  13  for  the  range 


of  ka  between  1.8  and  3.0.  There  are  two  well-defined  resonant  maxima 


shown  in  that  range  for  a positon  r/a  = 0.  934,  A = 45  in  the  first  basin 
(A)  whereas  for  the  center  of  the  second  basin  (B)  (Fig.  14)  only  one. 
rather  poorly  defined  maximum,  is  apparent  for  the  same  wave  number 


The  effect  of  location  can  be  eliminated  when  investigating  resonance 
by  plotting  the  maximum  amplification,  regardless  of  location,  as  a function 


of  ka.  The  theoretical  response  curve  defined  in  that  way  is  presented 
in  Fig.  16  for  these  two  coupled  circular  basins.  The  eight  resonant 
modes  are  immediately  obvious.  To  see  where  these  modes  of  oscillation 


come  from,  the  maximum  response  for  a circular  harbor  with  a 10° 
entrance  gap  directly  connected  to  the  open-sea  and  a circular  harbor 
with  a 60°  opening  directly  connected  to  the  open-sea  has  been  deter- 
mined using  the  "arbitrary  harbor  theory"  presented  by  Lee  (1969). 

These  theoretical  results  are  presented  in  Fig.  17  and  show  an  increase 
in  the  amplification  at  resonance  as  the  harbor  entrance  width  decreases 
along  with  a corresponding  shift  of  the  maxima  to  smaller  wave  numbers. 
The  former  is  referred  to  as  the  "harbor  paradox"  by  Munk  & Miles  (1961). 
The  latter  effect  is  the  shift  of  resonance  toward  closed  basin  resonance 
as  the  entrance  width  decreases. 

One  effect  of  coupling  can  be  seen  clearly  by  comparing  Fig.  16 
and  Fig.  17.  At  first  glance,  Fig.  16  appears  to  be  simply  a combination 
of  the  two  response  curves  presented  in  Fig.  17,  i.  e.  , the  first,  third, 
fifth,  and  seventh  modes  in  Fig.  16  correspond  to  the  four  modes  shown 
by  solid  lines  in  Fig.  17  and  the  second,  fourth,  sixth,  and  eighth  modes 
of  Fig.  16  agree  qualitatively  with  the  four  modes  by  dashed  lines  in 
Fig.  17.  Thus,  it  appears  that  the  maximum  response  of  this  coupled- 
basin  system  is  a combination  of  the  response  of  the  individual  basins 
when  each  is  connected  directly  to  the  open-sea.  This  suggests  one 
problem  which  can  arise  in  an  actual  harbor  when  it  is  enlarged  by 
changing  its  configuration  from  a single  basin  to  a system  of  inter- 
connected slips.  The  original  response  of  the  harbor  may  be  changed 
to  a more  complicated  one  with  additional  maxima  for  the  same  range 
of  wave  period.  (However  it  is  possible  that  this  effect  may  not  be  too 
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serious  due  to  the  influence  of  viscous  dissipation  which  tends  to  reduce 
the  maximum  amplification  for  modes  of  oscillation  whose  peaks  have 
small  wave  number  bandwidths.) 

The  shape  of  the  water  surface  for  the  eight  resonant  modes  of 
oscillation  shown  in  Fig.  16  are  presented  in  Figs.  18  through  25.  These 
may  be  compared  to  similar  distributions  for  the  individual  basins  (10° 
and  60°  openings)  which  have  been  presented  previously  by  Lee  (1969). 

Fig.  18  shows  the  wave  amplitude  distribution  inside  the  harbor 
for  the  first  resonant  mode  of  the  coupled-basins  system  (ka  = 0.  26). 

The  wave  amplitude  has  been  normalized  with  respect  to  that  at  location 
P (in  the  inner  basin)  which  is  a position  at  which  the  amplitude  is  very 
close  to  a maximum  (the  maximum  is  at  the  back-wall  on  the  diameter 
of  symmetry).  The  wave  amplitude  shown  in  Fig.  18  is  relatively  uniform 
within  each  basin;  however,  the  average  wave  amplitude  in  the  inner  basin 
is  about  twice  that  in  the  outer  basin.  In  addition,  either  positive  or  negative 
water  surface  displacements  occur  simultaneously  in  both  the  inner  and 
outer  basins.  Therefore,  this  mode  is  usually  referred  to  as  the 
"pumping  mode". 

The  wave  amplitude  distribution  inside  the  harbor  for  the  second 
mode  of  oscillation  (ka  = 0.  72)  is  presented  in  Fig.  19.  The  maximum 
wave  amplitude  is  located  in  the  outer  basin;  the  wave  amplitude  shown 
has  been  normalized  with  respect  to  that  at  location  P (which  is  very  close 
to  the  position  of  the  maximum  wave  amplitude).  Fig.  19  shows  that 
positive  water  surface  displacements  exist  in  the  outer  basin  while 


3.150);  Fig.  23  Shape  of  Mode  6 (ka  = 3.  450); 
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negative  water  surface  displacements  exist  in  the  inner  basin  with  a 
nodal  line  (line  of  zero  wave  amplitude)  located  near  the  common 
boundary  between  the  outer  and  inner  basins.  Considering  the  out-of- 
phase  motion  in  the  two  basins,  this  mode  of  oscillation  may  be  called 
the  "sloshing  mode",  i.  e.  , sloshing  between  basins.  The  increased 
complexity  of  water  surface  oscillations  with  increasing  wave  number 
is  seen  by  comparing  Figs.  18  and  19. 

Fig.  20  shows  the  wave  amplitude  distribution  inside  the  harbor 
for  value  of  ka  = 1.  94  which  corresponds  to  the  third  resonant  mode  of 
oscillation.  The  normalized  surface  amplitude  exhibits  two  nodal  lines: 
one  in  the  outer  basin  and  one  in  the  inner  basin.  Thus,  for  this  value 
of  ka  a sloshing  mode  exists  in  each  basin  with  the  water  surface  at  the 
inner  part  of  the  outer  basin  and  the  outer  part  of  the  inner  basin  moving 
in  phase.  In  addition,  for  this  mode  the  water  surface  displacement  in 
the  inner  basin  is  at  least  twice  the  wave  amplitudes  in  the  outer  basin. 

The  wave  amplitude  distribution  inside  the  harbor  for  the  fourth 
resonant  mode  (ka  = 2.34)  is  presented  in  Fig.  21.  The  position  of  the 
maximum  wave  amplitude  for  this  mode  is  in  the  outer  basin  near  the 
entrance  to  the  inner  basin.  There  are  three  nodal  lines  in  the  harbor: 
one  in  the  outer  basin,  one  near  the  common  boundary  between  the  two 
basins  (the  entrance  to  the  inner  basin),  and  another  near  the  center  of 
the  inner  basin.  Even  though  the  water  surface  displacements  in  the 
inner  basin  are  relatively  small  compared  with  those  in  the  outer  basin, 
the  water  surface  shape  for  this  value  of  ka  is  still  approximately  that 


of  a sloshing  mode  within  each  basin.  However,  this  mode  of  oscillation 
is  more  complicated  than  the  one  shown  in  Fig.  20  in  the  sense  that  an 
additional  nodal  line  exists  near  the  common  boundary  between  the  two 
basins. 

A wave  amplitude  distribution  for  the  next  resonant  mode  of 
oscillation  (ka  = 3.  15)  is  presented  in  Fig.  22.  For  this  mode  the 
maximum  wave  amplitude  occurs  in  the  inner  basin  at  the  position 
r = a,  9 = 90°  (and  its  symmetric  counterpart  at  9 = 270°).  There 
are  two  nodal  lines  in  each  basin  with  the  nodal  lines  somewhat  similar 
to  the  crossed  nodal  lines  of  the  corresponding  mode  of  oscillation  for 
a closed  circular  basin.  The  mode  shape  for  the  inner  basin  is  similar 
to  the  one  shown  in  Fig.  6.23  of  Lee  (1969),  while  the  mode  shape  for 
the  outer  basin  is  similar  to  that  presented  in  Fig.  6.24  of  Lee  (1969) 
except  near  the  entrance  to  the  inner  basin.  As  can  be  seen  from  Fig. 

22  water  surface  displacements  in  the  inner  basin  are  significantly 
larger  than  those  of  the  outer  basin. 

The  mode  shape  for  ka  = 3.45  is  presented  in  Fig.  23.  This 
value  of  ka  corresponds  to  the  sixth  resonant  mode  of  oscillation  shown 
in  Fig.  16.  For  this  mode,  the  maximum  wave  amplitude  occurs  in  the 
outer  basin  at  the  position:  r = a,  9 = 100^  (and  its  symmetric 

counterpart  9 = 260°).  There  are  two  nodal  lines  in  the  outer  basin, 
another  one  near  the  entrance  to  the  inner  basin,  and  two  nodal  lines  in 
the  inner  basin.  Comparing  Figs.  22  and  23  it  is  seen  that  the  water 
surface  shapes  appear  similar;  however,  corresponding  regions  in  the 
basins  are  in  phase  in  Fig.  22  (ka  = 3.  15)  while  they  are  out  of  phase 
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The  wave  amplitude  distribution  inside  the  harbor  for  the  seventh 
resonant  mode  of  oscillations  (ka  = 3.86)  is  presented  in  Fig.  24.  For 
this  value  of  ka,  the  maximum  wave  amplitude  occurs  at  the  center  of  the 
inner  basin,  and  the  water  surface  displacements  in  the  inner  basin  are 
much  larger  than  those  in  the  outer  basin.  The  contour  lines  of  constant 
amplitude  near  the  center  of  the  inner  basin  (including  the  nodal  line) 
are  nearly  circular.  It  is  also  seen  from  Fig.  24  that  water  surface 
motions  in  the  regions  near  the  center  of  both  the  inner  and  outer  basins 
are  in  phase. 

Fig.  25  shows  the  wave  amplitude  distribution  inside  the  harbor 
for  ka  = 4.  05  (the  eighth  resonant  mode  of  oscillation).  The  position  of 
maximum  wave  amplitude  is  at  the  center  of  the  outer  harbor,  and  the 
shape  of  the  water  surface  in  the  outer  basin  is  similar  to  that  shown  in 
Fig.  6.  26  of  Lee  (1969)  for  a circular  harbor  with  a 60°  opening  coupled 

» 

directly  to  the  open-sea.  For  this  wave  number  water  surface  displace- 
ments in  the  inner  basin  are  small  compared  with  those  of  the  outer  basin 
and  in  contrast  to  the  distribution  shown  in  Fig.  24  the  water  surface 
oscillations  near  the  center  of  the  two  basins  are  out-of-phase. 

Thus,  although  it  is  seen  that  the  response  of  coupled  basins  and 

I 

the  amplitude  distributions  within  a harbor  can  become  quite  complicated 
as  individual  basins  are  connected,  it  may  be  possible  to  investigate  these 
facets  of  the  problem  by  viewing  the  basins  as  separate  units.  For  example, 
in  the  case  of  the  coupled -circular  harbor  just  discussed  some  aspects  of 
the  response  and  the  water -surface  displacements  could  be  constructed  from 
a knowledge  of  the  response  of  the  individual  basins. 
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4.  3 THE  RESPONSE  OF  VARIOUS  COUPLED- BASINS  SYSTEMS 


In  Section  4.  I,  it  has  been  shown  that  the  response  of  the  Long 


Beach  Harbor  Model  obtained  using  the  coupled-basins  theory  agrees 


well  with  the  theoretical  results  (single-basin  theory)  and  experimental 


data  obtained  by  Lee  (1969).  This  agreement  demonstrates  how  well 


the  present  method  can  be  applied  to  basins  of  arbitrary  shape  coupled 


together  and  excited  by  waves  from  the  open- sea.  In  Section  4.  2, 


experimental  data  for  coupled -circular  basins  have  also  provided 


further  evidence  in  support  of  the  theoretical  approach,  although  in  the 


hydraulic  model  the  effect  of  viscous  dissipation  does  significantly 


reduce  the  amplification  factor  at  resonance,  especially  for  peaks  where 


the  wave  number  bandwidth  is  small. 


In  this  section,  the  effect  on  the  overall  harbor  response  of 


coupling  basins  together  will  be  explored  analytically  for  various  basin 


arrangements.  For  simplicity  rectangular  basins  and  circular  basins 


of  different  radii  and/or  different  openings  are  used.  The  purpose  of 


this  portion  of  the  investigation  was  twofold:  to  show  the  applicability 


of  the  method  developed,  and  to  study  the  effect  of  coupling  on  harbor 


response.  For  this  reason  only  the  maximum  response  curve  for  each 


case  will  be  presented  and  discussed. 


4.  3.  1 Two  Coupled  Circular  Harbors 


In  Fig.  26,  the  maximum  response  curve  is  presented  for 


a circular  harbor  with  a 10  opening  connected  along  its  diameter  of 


symmetry  to  another  circular  harbor  of  the  same  size  also  with  a 10C 


I 
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opening.  As  before,  the  ordinate  represents  the  maximum  amplification 
factor  defined  as  the  maximum  wave  amplitude  inside  the  harbor  (regard- 
less of  the  location)  divided  by  the  standing  wave  amplitude  at  the  harbor 
entrance  when  the  entrance  is  closed.  The  abscissa  is  the  wave  number 


parameter,  ka,  where  a is  the  radius  of  the  outer  harbor  (the  same  as 
the  inner  harbor  for  this  case).  From  Fig.  26  it  is  seen  that  there  are 
eight  maxima  in  the  range  of  ka  presented  corresponding  to  eight 


resonant  modes  of  oscillation.  The  values  of  ka  for  these  resonant 


position  at  which  the  wave  amplitude  reaches  the  maximum  is  located  in 


Basin  B (inner  basin)  for  the  1st,  3rd,  5th  and  7th  resonant  modes  of 


oscillation,  and  for  the  2nd,  4th,  and  6th  and  8th  resonant  modes  the 


maximum  is  in  Basin  A (outer  basin).  The  location  of  each  maximum  as 


being  in  either  Basin  A or  B is  indicated  in  the  figure  along  the  abscissa. 
The  effect  of  coupling  on  the  harbor  response  can  be  inferred  by  comparing 


Fig.  26  with  the  response  curve  for  a single  circular  harbor  of  10 
opening  connected  to  the  open-sea,  presented  in  Fig.  17.  The  response 
of  these  coupled  basins  appears  to  be  composed  of  the  resonant  modes 
of  the  individual  circular  harbors  with  some  differences  in  both  the 


resonant  wave  number  and  the  maximum  amplification.  Thus,  Fig.  26 


demonstrates  that  the  number  of  resonant  modes  of  oscillation  for 


a coupled  basin  is  larger  than  the  number  of  resonant  modes  for  either 


basin  alone,  and  for  two  identical  basins  connected  the  number  of  modes 


is  increased  by  about  a factor  of  two  for  a given  wave  number  range 


the  harbor  is  demonstrated  from  a comparison  of  Fig.  26  and  Fig.  16; 
in  the  latter  the  outer  harbor  has  a 60°  opening.  It  is  seen  that  the 
maximum  amplifi cation  factor  for  most  resonant  modes  for  the  example 
shown  in  Fig.  16  is  less  than  that  for  the  corresponding  modes  shown 
in  Fig.  26;  however,  in  contrast  the  wave-number-bandwidth  for  the 
resonant  modes  in  Fig.  16  is  generally  greater  than  those  of  Fig.  26. 

These  results  again  show  that  as  the  entrance  width  decreases  the 
maximum  wave  amplifi  cation  within  the  harbor  increases  while  the  wave 

j 

number  bandwidth  near  resonance  decreases.  However,  the  reader  is 
reminded  that  for  the  prototype,  viscous  effects,  which  are  not  con- 
sidered in  the  present  theory,  may  be  very  important  in  limiting  the 
maximum  amplification  at  resonance  especially  for  resonant  peaks  with 
narrow  bandwidths.  CThis  has  been  demonstrated  by  the  experimented 
results  presented  in  Figs.  12  through  15,  and  as  mentioned  will  be 
discussed  in  Section  4.  4.) 

The  response  of  two  circular  harbors  with  the  entrances 
oriented  at  90°  to  one  another  was  investigated.  The  basins  have  the 
same  diameter  and  each  entrance  has  an  included  angle  of  10°.  The 
maximum  response  curve  for  this  harbor  system  is  presented  in  Fig.  27 
where  the  ordinate  and  abscissa  are  defined  earlier.  There  are  eight 
maxima  (four  in  each  basin)  in  the  response  curve  for  the  range  of  ka 
presented  corresponding  to  eight  resonant  modes  of  oscillation.  The 
values  of  ka  for  these  resonant  modes  are:  ka  = 0.21,  0.62,  1.98, 

2.  10,  3.  12,  3.  35,  3.  85,  and  3.  90.  It  is  noted  that  the  number  of 
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resonant  modes  of  oscillation  and  the  associated  wave  numbers  for  the 
arrangement  shown  in  Fig.  27  are  similar  to  the  response  curve  pre- 
sented in  Fig.  26  suggesting  that  for  these  two  harbors  the  coupling 
pattern  does  not  alter  the  harbor  response  significantly.  This  is 
reasonable  since  this  constant-depth  circular  harbor  system  has  the 
same  entrance  width  and  only  one  characteristic  dimension  (the  dia- 
meter of  the  harbor).  Thus,  the  position  of  the  boundary  opening  pro- 
bably will  not  affect  the  response  of  the  harbor  very  much  for  this  range 
of  wave  lengths. 

From  Fig.  27  it  is  seen  that  the  first  two  modes  of  oscillation 
are  almost  identical  to  the  first  two  resonant  modes  shown  in  Fig.  26. 
This  is  reasonable  since  the  basin  dimensions  are  the  same  for  the  two 
cases  but  only  the  arrangement  is  different.  For  a small  ratio  of  radius 
to  wave  length  the  effec'  on  the  harbor  response  of  the  arrangement 
should  not  be  too  important.  The  wave  amplitude  inside  Basin  B for 
the  first  mode  (the  pumping  mode)  is  about  twice  the  wave  amplitude 
in  Basin  A,  similar  to  the  case  shown  in  Fig.  26.  The  theoretical 
amplification  at  resonance  for  the  other  modes  present  in  this  wave 
number  range  does  vary  somewhat  with  basin  arrangement  but  not 
significantly. 

The  maximum  response  curve  for  a circular  harbor  with  a 10° 
opening  coupled  to  another  circular  harbor  with  one -half  the  diameter 
is  presented  in  Fig.  28.  The  larger  harbor,  which  is  connected  to  the 
open-sea,  has  a main  entrance  of  10°  and  a 10°  opening  in  the  backwall, 
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thus,  the  smaller  basin  has  a 20°  opening.  In  the  abscissa  of  Fig.  28 
the  characteristic  dimension  "a"  is  the  radius  of  the  larger  basin 
(Basin  A). 

Fig.  28  shows  six  maxima  corresponding  to  the  six  resonant 
modes  of  oscillation  for  the  range  of  ka  presented;  the  values  of  ka 
for  these  are:  ka  = 0.  30,  0.  92,  2.  12,  3.  24,  3.  84  and  4.  06.  The 

position  of  the  maximum  amplitude  for  the  first  mode  (the  pumping  mode) 
is  at  the  center  of  the  backwall  of  Basin  B as  it  was  with  the  other  cases 
presented  in  the  previous  sub- sections.  It  is  interesting  to  note  that  the 
position  at  which  the  maximum  amplification  is  reached  is  in  Basin  A 
(the  larger  basin)  for  the  third,  fourth,  and  fifth  resonant  modes,  and 
for  the  first,  second,  and  sixth  resonant  modes  the  maximum  ampli- 
fication is  in  Basin  B.  The  wave  numbers  of  the  modes  at  ka  = 0.  30, 

2.  12,  3.24,  and  3.84  are  similar  to  those  of  the  four  modes  presented 
in  Fig.  17  for  a circular  harbor  with  a 10°  opening.  If  the  radius  of  the 
smaller  basin  (Basin  B)  is  used  as  the  normalizing  dimension  the  wave 
number  parameters  for  the  two  modes  at  ka  = 0.  92  and  4.  06  become 
ka  = 0.46  and  2.03  which  are  slightly  greater  than  the  values  for  the 
"pumping"  and  "sloshing"  modes  respectively  for  the  single  10°  harbor. 
This  is  in  the  right  direction  since  the  inner  harbor  has  a 20°  opening. 
Thus  the  results  presented  in  Fig.  28  further  demonstrate  the  fact  that, 
for  these  shapes,  the  resonant  modes  of  oscillation  of  coupled-basins 
are  combinations  of  the  resonant  modes  of  the  individual  basins  for  the 
range  of  ka  considered.  Therefore,  connecting  one  harbor  to  another 
does  not  necessarily  improve  the  amplification  characteristics  of  the 
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harbor  system;  however,  it  is  possible  that  the  wave  amplitude  ampli- 
fication for  some  particular  modes  may  be  reduced  if  the  proper 
geometry  is  found  and  if  viscous  effects  are  considered. 

The  maximum  response  curve  for  a circular  harbor  with  a 60° 
opening  coupled  to  a circular  harbor  ox  one-half  that  size  with  a 20° 
opening  is  presented  in  Fig.  29.  Similar  to  the  previous  example,  the 
characteristic  dimension  "a"  which  is  used  is  the  radius  of  the  larger 
basin  (Basin  A).  From  Fig.  29  it  is  seen  that  there  are  five  resonant 
modes  of  oscillation  for  the  range  of  ka  considered;  the  values  of  ka 
for  these  modes  are:  ka  = 0.41,  0.97,  2.32,  3.39  and  3.88.  For 
the  first,  second,  and  fifth  modes  (ka  = 0.41,  0.97,  and  3.88)  the 
maximum  wave  amplitude  occurs  in  the  smaller  basin  (Basin  B)  while 
for  the  modes  at  ka  = 2.  32,  3.  39  (third  and  fourth)  the  maximum 
response  is  in  the  larger  basin  (Basin  A).  By  comparing  Fig.  29  to 
Fig.  28  it  is  seen  that  except  for  the  mode  at  ka  = 3.  88  the  maximum 
amplification  factor  at  resonance  is  smaller  and  the  response  near 
resonance  has  a larger  wave  number  bandwidth  for  the  example  where 
the  outer  harbor  has  a 60°  opening  than  for  the  one  with  a 10°  opening. 
Thus,  the  effect  of  entrance  width  on  the  response  of  harbors  to  incident 
waves  is  further  demonstrated. 

The  results  presented  in  this  section  have  demonstrated  that  the 
modes  of  oscillation  which  exist  in  coupled-circular-basins  are  closely 
related  to  the  modes  of  the  individual  basins,  and  that  resonance  in  one 
basin  will  affect  the  oscillations  in  the  other.  Moreover,  for  a particular 
range  of  ka,  the  modes  of  oscillation  in  coupled-basins  appear  to  be  a 
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combination  of  the  possible  modes  of  oscillation  of  the  individual  basins. 

This  is  in  agreement  with  the  earlier  work  of  Raichlen  and  Ippen  (1965) 
in  which  it  was  found  that  a rectangular  harbor  connected  to  a highly 
reflective  and  larger  rectangular  wave  basin  had  many  more  modes  of 
oscillation  compared  to  the  response  of  the  same  harbor  connected 
directly  to  the  open-sea. 

4.  3.  2 Rectangular  and  Circular  Coupled  Harbors 

In  the  previous  subsection  examples  of  coupled  circular 
harbors  have  been  considered;  in  this  subsection  the  response  of  circular 
harbors  connected  to  rectangular  entrance  channels  will  be  discussed. 

I 

The  entrance  channel  is  a common  feature  of  prototype  harbors  and  the 
influence  of  the  channel  on  the  harbor  response  is  therefore  an  important 
feature  in  the  design  of  harbors. 

The  harbor  model  which  was  chosen  for  the  theoretical  investi- 
gation consists  of  a circular  basin  of  constant  depth  with  a radius  of 
0.  75  ft  and  an  entrance  gap  of  10°  included  angle  (an  entrance  width  of 
0.  131  ft).  Connected  to  this  entrance  is  a channel  0.  198  ft  wide  with  a 
variable  length.  The  other  end  of  the  channel  is  fully  open  and  communi- 
cates directly  with  the  open-sea.  There  is  a difference  in  the  width  of 
the  channel  compared  to  the  entrance  width;  the  channel  width  was  chosen 
to  be  the  same  as  the  rectangular  channel  investigated  previously  and 
reported  by  Lee  (1969)  and  the  dimension  of  the  circular  harbor  were 
the  same  as  that  discussed  previously  herein. 
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The  response  curves  for  the  maximum  response  anywhere  in 
the  harbor  are  presented  in  Figs.  30  through  33  for  four  different  lengths 
of  the  entrance  channel.  These  response  curves  were  obtained  by  using 
the  theory  presented  in  Chapter  2 by  dividing  the  entire  domain  into 
3 regions:  the  open-sea  region,  the  entrance  channel  region  (Basin  A), 

and  the  circular  harbor  region  (Basin  B).  In  all  of  these  figures  the 
ordinate  is  the  maximum  amplification  factor  in  the  harbor  regardless 
of  location  and  the  abscissa  is  the  product  of  the  wave  number  and  the 
radius  of  the  circular  basin.  Indicated  near  the  abscissa  in  these  figures 
is  the  particular  basin  where  the  maximum  occurs,  e.  g.  , "A"  indicates 

the  maximum  occurs  in  the  entrance  channel  and  "B"  denotes  a maximum 
in  the  circular  harbor.  Considering  the  shortest  entrance  channel  to  be 
of  unit  length, the  other  three  lengths  studied  were  two,  three,  and  four 
times  that  length. 

The  response  with  the  shortest  entrance  channel,  a length  approxi- 
mately one-third  the  diameter,  is  presented  in  Fig.  30.  Five  resonant 
modes  of  oscillation  are  evident  for  the  range  of  ka  investigated:  three 

occur  in  the  circular  basin,  one  in  the  entrance  channel,  and  one  mode  has 
maximum  amplitudes  in  both  basins.  Referring  to  Fig.  17,  it  is  noted 
that  a circular  harbor  with  a 10°  opening  has  four  modes  of  oscillation 
over  a range  of  ka  up  to  about  4.  0.  The  values  of  the  wave  number  para- 
meter for  the  first,  second,  fourth  and  fifth  modes  of  Fig.  30  are  about 
the  same  as  those  for  the  four  modes  of  oscillation  of  the  10°-opening- 
harbor  shown  in  Fig.  17.  The  amplification  at  resonance  is  different, 
in  fact  the  amplification  for  the  pumping  mode  is  increased  by  nearly 
45%  by  adding  this  fairly  short  channel. 
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Fig.  32  Maximum  Response  Curve  for  Circular  Harbor  Coupled  to 
Rectangular  Entrance  Channel  (l/a.  = 2.05) 
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Fig.  33  Maximum  Response  Curve  for  Circular  Harbor  Coupled  to 
Rectangular  Entrance  Channel  (-t/a  = 2.73) 
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When  the  length  of  the  entrance  channel  is  doubled,  as  is  shown 
in  Fig.  31,  an  additional  maximum  is  introduced  in  the  response  compared 
to  Fig.  30  for  the  same  wave  number  range.  The  four  modes  of  the  cir- 
cular harbor  become  quite  distinct  and  comparable  in  wave  number  to 
those  shown  in  Fig.  17  for  the  same  circular  harbor  connected  directly 
to  the  open- sea.  The  two  maxima  shown  in  Fig.  31  at  ka  = 1.  58  and 
3.  67  correspond  to  maximum  amplitudes  in  the  entrance  channel.  By 
multiplying  these  values  by  the  ratio  of  the  channel  length  to  the  harbor 
radius,  -t/a,  the  wave  number  parameter  may  be  expressed  in  terms  of 
the  channel  length  as:  k-C  = 2.  14  and  4.  83.  The  value  of  this  parameter 

for  the  first  two  modes  of  a rectangular  harbor  (with  the  same  aspect 
ratio)  connected  directly  to  the  open-sea  is:  kf-  = 1.32  and  4.2.  There- 

fore, when  the  rectangular  channel  is  connected  between  the  open-sea 
and  the  circular  basin  the  tendency  is  for  the  mode  of  oscillation  to  shift 
toward  the  closed  rectangular  basin  mode  (k-t  = it  and  2tt).  Nevertheless, 
these  are  identifiable  as  modes  of  resonance  for  the  entrance  channel. 

Fig.  32  shows  the  maximum  response  curve  when  the  length  of 
the  entrance  channel  is  increased  by  50%  compared  to  the  case  shown 
in  Fig.  31.  For  this  configuration  the  four  modes  of  the  circular  basin 
can  still  be  recognized  and  are  located  at  approximately  the  same  values 
of  ka  as  before,  but  the  wave  numbers  for  the  modes  of  oscillation  in 
the  entrance  channel  have  changed  and  the  number  of  modes  in  the  channel 
have  now  increased  to  three.  Thus,  there  are  now  seven  modes  of 
oscillation  for  the  range  of  the  wave  numbers  shown,  i.  e.  , 0 < ka  < 4.  0. 
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basin;  however,  the  ratios  of  channel  length  to  wave  lengths  for  reson- 
ance are  comparable  to  the  closed  basin  and  open  harbor  resonant  modes 
one  would  expect  for  the  rectangular  basin  alone.  Another  trend  which  is 
evident  is  that,  for  the  four  cases  chosen,  the  amplification  at  resonance 
for  the  pumping  mode  increases  with  increasing  channel  length;  for  ratios 
of  channel  length  to  harbor  radius  of  0.  68,  1.  35,  2.  05,  and  2.  74  the 
corresponding  values  of  the  amplification  of  the  pumping  mode  are  8.  5, 

9.  1,  13.95,  and  17.0.  This  is  compared  to  a value  for  the  circular 
harbor  alone  of  5.9.  Thus,  again  the  response  of  a more  complicated 
harbor  appears  qualitatively  to  be  a combination  of  the  response  of  the 
individual  basins  with  certain  effects  relating  to  the  influence  of  one 
basin  type  on  the  other. 

Various  investigators  have  proposed  the  use  of  side  channel 
resonators  to  eliminate  or  at  least  reduce  the  effect  of  resonance  on 


harbors  (see  Valembois  (1953),  James  (1968)).  These  resonators 
are  rectangular  channels  of  various  lengths  which  are  connected  per- 
pendicular to  the  main  entrance  channel.  Since  the  effect  of  the  length 
of  an  entrance  channel  on  the  response  of  a circular  harbor  had  been 
investigated,  it  was  considered  ioglv.c.1  to  extend  this  to  a minimal 
investigation  of  the  effect  of  resonators  of  two  different  lengths  on  the 
response  of  the  harbor  shown  in  Fig.  33.  The  channels  were  located 
approximately  one -third  of  the  length  of  the  entrance  channel  from  the 
open- sea  with  a width  equal  to  the  width  of  the  entrance  channel.  The 
cases  chosen  were  for  lengths  of  the  resonators  equal  to  one -half  and 
one -quarter  of  the  length  of  the  main  channel.  The  maximum  response 
curves  for  these  two  harbor  systems  are  presented  in  Figs.  34  and  35. 
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Fig.  34  Maximum  Response  Curve  for  Circular  Harbor  Coupled  to 

Rectangular  Entrance  Channel  with  Side  Chambers  (c/’t.  = £) 
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Fig.  34  shows  the  response  curve  presented  with  the  ordinate  and 
abscissa  as  previously  defined  for  the  case  with  the  resonators  one-half 
the  length  of  the  main  entrance  channel;  the  arrangement  is  shown  in  the 
inset  in  this  figure.  For  this  harbor  system  the  wave  number  of  the 
pumping  mode  has  changed  only  slightly  while  its  maximum  response  has 
been  reduced  from  about  seventeen  to  ten  by  the  addition  of  the  resonator. 
The  next  mode  has  been  shifted  to  a smaller  wave  number  and  attenuated 
by  a factor  of  nearly  two  by  the  resonator.  It  is  noted  that  two  curves  are 
shown  for  this  peak:  one  where  resonance  is  in  the  entrance  channel  (the 
solid  curve)  and  a curve  for  which  resonance  is  in  the  side  channel  or 
resonator  (the  curve  composed  of  long  dashes).  For  this  case,  resonance 
in  the  side-channel  dominates;  however,  with  respect  to  the  overall 
response  of  the  harbor  this  is  not  considered  important.  The  next  four 
modes  are  not  modified  appreciably  except  for  some  shift  in  the  resonant 
wave  numbers.  The  seventh  mode  which  is  associated  with  resonance  in 
the  circular  harbor  is  affected  by  the  side  channel  resonators  by  both  an 
increase  in  amplification  compared  to  the  case  shown  in  Fig.  33  and  by  a 
decrease  in  the  bandwidth  of  the  peak.  Although  theoretically  this  appears 
to  result  in  a more  responsive  harbor,  the  reduction  in  the  bandwidth  of 
the  mode  makes  it  more  susceptible  to  viscous  effects  and  hence  the  actual 
response  may  be  less  than  for  the  corresponding  mode  shown  in  Fig.  33 
(see  Section  4.  4). 

The  response  curve  which  is  obtained  when  the  length  of  the  side 
channel  resonator  is  halved  is  presented  in  Fig.  35.  The  peaks  which  are 
labeled  as  B (occurring  in  the  circular  basin)  are  similar  in  all  three 
cases  (Figs.  33,  34,  and  35).  There  is  a slight  change  in  the  maximum 


Fig.  36  Response  Curve  for  a Circular -Segment  Harbor  with 
Entrance  Channel;  at  Center  of  Backwall  (r  = a,  9 


into  two  regions:  the  rectangular  entrance  channel  and  the  circular 

sector.  The  value  of  ka  for  the  first  two  modes  of  oscillation  are: 
ka  = 0.289  and  ka  = 3.  744.  The  first  resonant  mode  of  oscillation  is 
a pumping  mode  where  the  maximum  amplification  (an  amplification 
factor  of  7.  48)  occurs  at  the  center  of  the  backwall.  This  value  of 
amplifi cation  factor  is  very  much  smaller  than  the  value  reported  by 
Carrier,  Shaw,  and  Miyata  (1971)  in  which  they  have  reported  an 
amplification  factor  of  about  22.  For  the  second  resonant  mode  the 
amplification  factor  computed  using  the  present  theory  is  2.  75  which 
is  close  to  the  value  of  about  2.  9 reported  by  Carrier,  Shaw,  and 
Miyata  (1971). 

Shown  as  insets  in  Fig.  36  are  the  computational  details  of  the 
response  near  resonance  for  the  two  modes;  the  data  from  the  coupled- 
basins  theory  are  indicated  by  solid  circles.  It  is  seen  that  the  maxi- 
mum amplification  was  effectively  determined  for  each  mode  of  oscillation. 
Thus,  it  appears  that  the  difference  in  the  amplification  factor  for  the 
pumping  mode  between  the  present  theory  and  the  theory  of  Carrier, 

Shaw,  and  Miyata  can  be  attributed  to  differences  between  the  methods 
rather  than  incomplete  computations  near  resonance. 

In  addition  to  the  theoretical  results,  some  experiments  were 
conducted  in  the  laboratory  for  comparison.  Due  to  certain  experi- 
mental limitations  and  to  cover  the  range  of  ka  for  the  first  mode  of 
oscillation  it  was  necessary  to  use  a small  harbor  model.  Referring  to 
the  sketch  shown  in  Fig.  36,  the  dimensions  of  the  harbor  model  used 
for  the  experiment  were:  9q  = 70°,  a = 3 in.  , £ = 1.  5 in.  , b = 0.  18 

in.  Experimental  data  are  shown  as  open  circles  in  Fig.  36.  It  is  seen 
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that  away  from  resonance  the  theory  agrees  reasonably  well  with  the 
experiments.  However,  at  resonance  the  maximum  amplification 
measured  was  approximately  40%  of  the  maximum  predicted  by  the 
linear  inviscid  theory;  this  reduction  was  due  undoubtedly  to  viscous 
effects  which  were  quite  important  in  these  experiments  because  of  the 
small  width  to  depth  ratio  of  the  entrance  channel. 

The  amplitude  distributions  determined  theoretically  for  the  first 
two  mode  of  oscillation  are  presented  in  Figs.  37  and  38.  Fig.  37  shows 
the  water  surface  elevation  for  the  pumping  mode  (ka  = 0.  289)  increasing 
radially  with  distance  from  the  entrance  with  the  wave  amplitudes  inside 
the  circular -soctor  region  nearly  uniform  in  the  9-direction  at  a given 
radius.  The  variation  of  the  water  surface  elevation  inside  the  harbor 
for  the  second  mode  of  oscillation  (ka  = 3.  744)  is  shown  in  Fig.  38. 

For  this  mode  of  oscillation  the  maximum  wave  amplitude  occurs  within 
the  entrance  channel  relatively  near  the  entrance  to  the  circular  segment 
region  with  a nodal  line  present  inside  the  circular  sector  basin.  Thus, 
negative  water  surface  displacements  are  produced  in  the  region  near  the 
back  wall. 

The  results  presented  in  this  chapter  have  shown  the  applicability 
of  the  present  theory  in  analyzing  the  coupled  basins  problem.  For  a 
complicated  harbor  geometry  it  may  not  be  possible  or  economical  to 
analyze  the  response  using  the  method  developed  by  Lee  (1969);  however, 
problems  of  computer  storage  size  and  economy  of  computation  may  be 
overcome  using  the  present  theory  where  the  harbor  region  is  divided 
into  several  basins.  The  boundary  segment  size  in  each  basin  must 
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still  be  small  compared  to  the  wave  length;  the  criterion  suggested  by 
Lee  (1969)  that  the  maximum  boundary  segment  size  be  smaller  than 
one-tenth  of  the  smallest  wave  length  studied  also  applies  to  the  present 
theory. 

4.  4 THE  EFFECT  OF  ENERGY  DISSIPATION  ON  HARBOR  RESONANCE 

In  this  section  the  effect  of  energy  dissipation  on  the  response  of 
a harbor  will  be  discussed  qualitatively  with  attention  given  to  the  differ- 
ences between  hydraulic  models  and  the  corresponding  prototype  harbors 
near  resonance. 

There  are  three  major  effects  of  viscous  dissipation  on  the  response 
of  a dynamic  system;  to  describe  these  it  is  useful  to  refer  to  the  simple 
example  of  the  forced  oscillation  of  a single-degree-of-freedom  oscillator 
(such  as  a spring-mass -dashpot  system)  described  by  the  following 
equation  of  motion: 

mx  + cx  + k „x  = F sin  u)t  (36) 

“ o 

where  m is  the  mass  of  the  oscillating  body,  c is  a coefficient  of  damping 
for  the  system,  k.  is  the  spring  constant,  x is  the  displacement,  and  F 

n*  0 

is  the  amplitude  of  an  applied  force  with  circular  frequency  uu.  Eq.  36 
may  be  rewritten  as: 

x + 2Qw  x + uuax  = X__,uj2  sin  uut  (37) 

n n ST  n ' 

where  C = c^cc>  cc  *s  defined  as  a critical  damping  coefficient  equal  to 
2muun  where  uu2  = k#An,  and  Xg,p  is  the  static  displacement  of  the  spring 
system  under  the  applied  force  Fq  and  is  equal  to  Fo/k*.  The  dynamic 


response  of  this  system  is  then  described  by: 


[‘-£)7*L«f] 


In  Eq.  38  M is  the  ratio  of  the  maximum  excursion  of  the  oscillating  mass 

at  the  given  frequency  to  the  static  movement  of  the  mass  caused  by  F , i. 

o 

Xg,p.  Eq.  38  is  the  well  known  expression  for  the  system  amplification 
which  indicates  that  the  amplifi cation  goes  to  infinity  for  zero  damping 
when  the  forcing  frequency  equals  the  natural  frequency,  i.  e.  , uo/u)n  = 1. 
The  phase  angle  between  the  forcing  function  and  the  system  output  is 


given  by: 


2£— 

ou 

n 


tan  <b  = 


-(£)' 


From  Eq.  38  it  is  seen  that  the  amplitude  at  resonance  can  be  expressed 
as  a first  approximation  for  small  damping  as: 

Mr  = 1/2C  , (40) 

where  the  subscript  R refers  to  resonance.  The  shift  of  the  resonant 
frequency  caused  by  damping  can  be  shown  to  be: 

(tu/u)  ) = J\ - ~2Qk  . (41) 


The  peakedness  of  a response  curve  is  usually  described  by  the 
frequency  bandwidth  of  the  half-power  point,  i.  e.  , the  point  at  which  the 
power  has  dropped  to  one -half  its  peak  value  or  the  amplitude  to  0.  707  of 


its  peak.  For  small  values  of  the  damping  factor,  C»  the  frequency 
limits  of  the  half- power  point  are  given  approximately  by: 


oi/ou  = J\  ± 2C  • (42) 

n 

which  gives  a bandwidth  for  the  half-power  point  of  A(u>/u;n)  2:  2Q. 

Eqs.  40,  41,  and  42  are  shown  in  Fig.  39  with  the  damping  factor 
as  the  abscissa  and  the  ordinate  as  the  amplification  factor  at  resonance, 
the  shift  of  the  resonant  frequency  due  to  damping,  and  the  frequency 
bandwidth  of  the  half-power  point.  One  obvious  feature  of  these  curves 
is  that  the  amplification  factor  and  the  bandwidth  of  the  half-power  point 
are  affected  much  more  by  damping  than  is  the  resonant  frequency.  For 
the  range  of  damping  factor  shown  (0.  05  < C,  < 0.  3)  the  amplification 
factor  at  resonance  and  the  half-power  bandwidth  change  by  a factor  of 
six  while  the  resonant  frequency  shifts  by  only  about  4%  to  lower  fre- 
quencies. Therefore,  for  the  single-degree-of-freedom  oscillator  the 
major  effect  of  increased  damping  on  the  resonant  response  are  to 
decrease  the  amplification  at  resonance  and  increase  the  frequency  band- 
width at  the  half-power  point  while  maintaining  approximately  the  same 
resonant  frequency  of  the  system. 

Miles  and  Munk  (1961)  and  Ippen  and  Goda  (1963)  have  shown  that 
the  dynamics  of  a harbor  near  resonance,  at  least  for  the  lowest  mode  of 
oscillation,  are  similar  to  the  single-degree-of-freedom  oscillator. 
Therefore,  in  qualitative  sense  one  would  expect  similar  effects  of 
damping  on  the  harbor  response  with  respect  to  the  amplification  at 
resonance,  the  frequency  bandwidth  of  the  half-power  point  and  the  shift 
in  the  resonant  frequency. 
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These  effects  would  apply  in  the  same  way  to  both  a hydraulic 
model  and  the  corresponding  prototype  harbor;  however,  the  relative 


energy  dissipation  (or  the  damping  factor)  in  the  two  cases  may  be 
quite  different,  i.  e.  , scale  effects  may  be  important.  In  both  model 
and  prototype  generally  there  are  at  least  five  regions  where  energy 
dissipation  can  be  important:  at  the  harbor  entrance,  at  the  internal 
boundaries  of  the  harbor,  along  the  bottom,  at  internal  structures  such 
as  piers,  moles  etc.  , and  internal  viscous  dissipation  in  the  wave 
system.  These  forms  of  dissipation  will  be  discussed  with  respect  to 
the  model  and  the  prototype,  and  some  attempt  will  be  made  to  indicate 
the  importance  of  scale  effects,  i.  e.  , the  reliability  of  scaling  up 
response  characteristics  of  the  model  to  the  prototype. 

The  dissipation  of  energy  at  the  harbor  entrance  can  be  divided 
into  two  parts:  the  energy  loss  associated  with  separation  and  that  due 
to  boundary  friction.  For  the  lower  modes  of  oscillation  where  wave 
periods  and  entrance  velocities  are  relatively  large  the  former  would 
probably  be  more  important.  If  the  effect  of  separation  is  important  in 
defining  the  entrance  losses  probably  it  is  equally  important  in  both  model 
and  prototype  and  scale  effects  for  this  type  of  dissipation  would  be 
relatively  unimportant.  Thus,  if  this  loss  could  be  determined  in  the 
model  it  could  be  directly  scaled  to  the  prototype.  The  effect  of  this 
loss  has  been  demonstrated  by  Ippen  and  Raichlen  (1962)  in  the  study  of 
the  response  of  two  highly  reflective  coupled  rectangular  basins.  In  that 
case  the  resonant  response  of  a small  rectangular  harbor  was  investigated 
with  and  without  a screen  (16  mesh/inch)  stretched  across  the  entrance. 
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geometrically  scaled  models  are  considered  the  value  of  the  ratio  of  a 
for  the  two  models  is  considered  to  be  an  indication  of  the  relative 
importance  of  dissipation  in  the  two  models.  The  subscript  r is  used 
to  denote  the  ratio  of  a quantity  in  the  smaller  model  to  the  corresponding 
quantity  in  the  larger  model.  For  the  same  relative  wave  number  in 
both  models,  i.  e.  , (kb)^  = (kh)f  = 1,  Eq.  44  reduces  to  the  following 
for  the  ratio  of  the  boundary  dissipation  moduli: 


a = t~ 
r b 


(45) 


Since  the  models  are  operated  as  Froude  models,  Tr 

L and  Eq.  45  becomes: 
r 


L^  and  b = 
r r 


Eq.  46  shows,  for  geometrically  similar  models,  the  smaller  model 
would  have  greater  relative  dissipation  than  the  larger  model.  However, 
if  the  models  are  not  geometrically  similar  then,  as  will  be  shown,  the 
complete  first  term  of  Eq.  44  must  be  used  in  making  a comparison. 

Lee  (1969)  presented  the  results  of  both  theory  and  experiment 
for  a model  of  the  East  and  West  Basins  of  Long  Beach  Harbor.  (These 
results  are  also  presented  in  Figs.  8 through  1 1 of  this  report  and 
compared  to  the  results  obtained  from  the  coupled -basins  theory.  ) 

Knapp  and  Vanoni  (1945)  used  ahydraulic  model  to  investigate  the 
alignment  of  a mole  to  protect  this  part  of  Long  Beach  Harbor  and 
some  selected  results  of  their  experiments  were  also  presented  and 


. . . . . . . . 
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discussed  by  Lee  (1969).  Lee's  experiments  were  conducted  in  a 
constant  depth  model  with  a horizontal  length  scale  of  1/4700  and  an 
average  vertical  length  scale  of  1/40.  The  model  of  Knapp  and  Vanoni 
more  realistically  modeled  depth  effects  (although  these  effects  are 

I 

very  small  for  this  harbor)  since  a horizontal  scale  of  1/480  and  a 
vertical  scale  of  1/240  was  used.  For  the  resonant  mode  of  oscillation 
which  corresponds  to  a period  of  about  6 min.  (ka  2:  3.  3,  where  "a" 
is  a typical  length  dimension  of  6768  ft  in  the  prototype)  it  was  found 
that  Lee's  model  gave  an  amplification  factor  nearly  2.  5 times  that 
observed  in  the  Knapp  and  Vanoni  model.  This  is  contrary  to  what  Eq. 

46  predicts  if  the  horizontal  scales  alone  are  considered;  for  these 
scales  Eq.  46  says  that  the  damping  modulus  for  the  smaller  harbor 
should  be  nearly  six  times  that  for  the  larger  model.  This  contradiction 
is  explained  by  the  fact  that  the  two  models  are  not  geometrically 
similar  and  the  complete  first  term  of  Eq.  44  must  be  used  in  making 
a comparison. 

The  characteristics  of  these  two  models  for  the  6 min.  mode  in 
the  prototype  are  shown  in  Table  1.  In  this  table,  the  parameter  "a" 
denotes  the  length  of  the  North  side  of  the  West  Basin  and  "b"  is  the 
maximum  width  of  that  basin.  Using  these  data  and  applying  Keulegan's 
analysis  for  the  sloshing  in  a rectangular  basin,  considering  only 
boundary  dissipation,  the  ratio  of  the  damping  modulus  a for  Lee's 
model  to  the  modulus  for  Knapp  and  Vanoni' s model  is  1/1240.  This 
shows  that  the  expected  dissipation  in  the  larger  model  would  be  much 


Lee 


Knapp  & Vanoni 


L 

1/4700 

1 /480 

rh 

L 

I* 

1/40 

1/240 

V 

a (ft) 

1.44 

14.  1 

b (ft) 

0.  95 

9.  27 

h (ft) 

1.  0 

0.  167 

k (1/ft) 

2.  26 

0.  23 

T (sec) 

0.  745 

11.  8 

greater  than  that  in  the  smaller  model  primarily  because  the  models  are 
not  geometrically  similar  and  the  vertical  scale  of  the  smaller  model 
is  more  exaggerated.  Obviously  the  magnitude  of  the  damping  modulus 
ratio  given  is  only  an  indication  of  the  relative  importance  of  dissipation 
at  the  two  scales  and  it  is  not  an  exact  measure  of  the  relative  dissipation 
in  the  two  models. 

A different  problem  arises  when  model  results  are  to  be  scaled 
up  to  prototype  systems,  since  the  boundary  layers  in  a model  are 
probably  laminar  whereas  in  the  prototype  turbulent  boundary  layers  are 
most  probable.  Since  Keulegan's  treatment  considers  laminar  damping 
only,  this  approach  is  not  applicable.  However,  to  determine  trends 
in  dissipation,  steady  flow  considerations  could  be  applied  to  this 
problem  for  the  lower  modes  of  oscillation,  once  the  wave  periods  of  these 
modes  are  usually  very  large.  Assuming  steady  flow,  the  boundary  drag 


A 


can  be  defined  in  terms  of  an  average  skin  friction  coefficient,  C^, 
which  takes  on  different  forms  for  laminar  and  turbulent  flow.  The 


power  dissipated  due  to  boundary  shear  can  be  expressed  as: 

Pd  = CfB2U3  (47) 

where  B is  a surface  area  and  U is  a characteristic  velocity  of  the 
system. 

If  the  model  were  truly  a Froude  model  the  power  dissipated 

7 /2 

would  scale  up  in  the  prototype  by  which  assumes  the  same  skin 

friction  coefficient  in  the  model  as  in  the  prototype,  i.  e.  , = 

r 

I Cj  = 1.  Since  the  model  is  assumed  to  operate  in  the  laminar 
m/  p 

flow  region  and  the  prototype  in  turbulent  flow,  this  ratio  will  not  be  unity 

and  the  magnitude  of  will  indicate  whether  proportionately  more 

r 

energy  is  dissipated  in  the  model  than  in  the  prototype,  e.  g.  , if  > l 

r 

the  model  would  dissipate  relatively  more  energy  than  the  prototype  and 

for  Cj  <1  the  opposite  is  true.  These  skin  friction  coefficients  can  be 
r 

expressed  as: 


for  laminar  flow: 


and  for  turbulent  flow: 


1.  33 


0.  074 


where  El  is  the  Reynolds  Number  based  on  a horizontal  length,  and  as 
before  the  subscripts  m and  p indicate  model  and  prototype  respectively. 
Hence,  for  a Froude  model  the  skin  friction  coefficient  ratio  becomes: 


• 'Mnimii  MRiMiwaNMMMMRNCHWEMB^ 
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wherein  L is  the  ratio  of  a horizontal  model  dimension  to  its  corres- 
r 

ponding  prototype  dimension.  For  the  sake  of  this  discussion  where 
only  trends  are  important  the  velocity  which  will  be  used  in  defining  the 
Reynolds  Number  is  the  maximum  water  particle  velocity  in  a wave  and 
the  wave  length  will  be  used  as  the  characteristic  length  dimension  cf 
the  system. 

It  is  conceivable  that  for  large  values  of  and  the  Reynolds 

Number,  would  be  less  than  unity  indicating  proportionately  more 

r 

dissipation  in  the  prototype  than  in  the  model.  However,  consider  the 
following  example:  the  length  scale  (h>r)  is  1/100,  the  prototype  depth 
is  30  ft,  and  the  wave  period  is  2.  8 min.  (the  wave  length  would  be 
5000  ft).  Taking  a wave  height  of  2 ft,  as  a first  approximation,  the 
maximum  water  particle  velocity  for  this  long  wave  would  be  1 ft/sec. 
The  Reynolds  Number  based  on  this  velocity  and  wave  length  would  be 

g 

about  5 x 10  . Therefore,  from  Eq.  48,  would  be  1.4,  and  for  this 

r 

case  proportionately  more  energy  would  be  dissipateu  on  the  bottom  in 
the  model  than  in  the  prototype. 

Two  aspects  of  energy  dissipation  in  harbors  have  not  been 
discussed  yet.  The  first  deals  with  the  damping  of  waves  inside  the 
harbor  due  to  the  permeability  of  the  bottom  and  the  second  deals  with 
energy  dissipation  at  the  boundaries  due  to  wave  run-up  and  breaking  on 
rough  embankments  and  structures  within  the  harbor. 
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I 

Reid  and  Kajuira  (1957)  have  treated  the  problem  of  viscous 

effects  associated  with  percolation  into  a permeable  sea  bed  of  infinite 

thickness  caused  by  temporal  and  spacial  pressure  distributions  at  the 

• Hv 

sea  bed.  Considering  exponential  damping  with  distance  (a  = aQe 
the  damping  modulus  given  by  them  is  at  worse: 


D = 0.  123  (49) 

vh  5 

and  less  for  values  of  h/LQ  greater  and  less  than  0.  13.  The  quantity 

-9 

f3  in  Eq.  49  is  a permeability  coefficient  equal  to  about  10~7  ft2  for 
ordinary  sand.  Thus,  for  a 30  ft  depth  Eq.  49  gives  a value  of  D of  a 
about  10  ^ ft“l  which  indicates  that  for  most  purposes  this  aspect  of 
dissipation  can  be  neglected  compared  to  other  forms  of  boundary 
dissipation. 

Energy  dissipation  around  the  harbor  periphery  due  to  wave  run 
up  and  breaking  on  beaches  and  revetted  structures  cannot  be  evaluated 
specifically  due  to  the  wide  range  of  structures  which  may  be  located 
at  the  boundaries  of  each  particular  harbor.  It  is  possible  that  with  the 
proper  design  of  internal  structures  the  energy  loss  in  a harbor  due  to 
this  type  of  dissipation  can  completely  overshadow  any  other  form  of 
energy  loss. 

Ippen  and  Goda  (1963)  analytically  investigated  the  effect  of 


energy  dissipators  on  the  response  of  a fully  open  narrow  rectangular 
harbor  connected  to  the  open-sea;  their  development  is  quite  instructive 
in  describing  the  effect  of  the  reflective  characteristics  of  the  structures 
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on  the  response  and  it  is  summarized  here.  They  determined  the 
response  of  the  harbor  as  a function  of  a reflection  coefficient  defined 
in  terms  of  the  standing  wave  amplitude  in  the  harbor  as: 


A = A0(l  +K) 


where  A is  the  standing  wave  amplitude,  Aq  is  the  amplitude  of  a pro- 
gressive wave  within  the  harbor,  and  K is  the  reflection  coefficient. 
Thus,  for  K = 1 (perfect  reflection)  A = 2Aq  and  for  K = 0 (zero 
reflection)  A = Aq.  When  the  response  is  determined  assuming  a two 
dimensional  oscillation  in  the  harbor  matched  at  the  entrance  to  the 
solution  for  the  open- sea,  the  following  simple  expression  can  be 
obtained  which  is  a first  approximation  to  the  response  of  the  narrow, 
fully-open  harbor  at  resonance: 


R - 1 _2_ 

R - kb + 6 * 1 - K 


(50) 


where  b is  the  width  of  the  harbor  and  is  equal  to  the  entrance  width  for 
the  fully-open  harbor  and  6 = (1  -K)/(l  +K).  As  an  example,  consider  a 
fully  open  harbor  where  the  ratio  of  the  width  of  the  harbor  to  its  length 
is  0.  02  and  the  normalized  resonant  wave  number  of  the  fundamental 
mode  is  kl  ^ ir/2.  From  Eq.  50  the  undamped  amplification  at  reson- 
ance would  be  63.  6.  For  boundary  reflection  coefficients  of  K = 0.  9, 

0.  7,  and  0.  5 the  corresponding  amplification  factors  at  resonance  are 
Rp  = 14.7,  5.2,  and  2.9  respectively.  Thus,  with  a general  internal 
reflection  coefficient  of  only  0.  9 the  harbor  response  would  be  decreased 
more  than  four  times  from  its  undamped  value. 


8 - 

In  the  case  of  a harbor  with  a more  complicated  shape,  at  the 
present  time  it  is  not  possible  to  evaluate  the  effect  of  this  type  of 
boundary  dissipation  on  the  response.  However,  the  example  just 
presented  indicates  that  this  effect  can  be  important  and  makes  it  very 
difficult  to  scale  up  the  results  of  a model  study  to  a prototype  scale 
with  any  reliability.  The  only  obvious  assurance  one  has  is  that  the 
response  would  be  less  than  an  inviscid  theory  predicts.  Nevertheless 
some  of  the  methods  describe  here  can  provide  a guide  to  an  engineering 
decision;  but  more  attention  must  be  given  to  this  aspect  of  the  problem 
in  future  research  to  lead  to  a more  generally  applicable  theory.  On 
the  other  hand,  as  mentioned  earlier,  the  inviscid  theory  developed  and 
presented  in  this  report  effectively  describes  the  important  wave 
periods  for  resonance  even  if  there  is  moderate  dissipation  in  the 
prototype;  this  in  itself  is  a useful  guide  for  engineering  design. 
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CHAPTER  5 
CONCLUSIONS 

The  following  major  conclusions  can  be  drawn  from  this  study: 

1.  The  theory  developed  herein  and  termed  the  coupled-basins 
theory  predicts  the  response  of  an  arbitrary  shape  harbor 
with  constant  depth  which  can  be  divided  into  several  inter- 
connected basins. 

2.  The  theory  agrees  well  with  experimental  results  obtained 
from  two  connected  circular  harbors  as  well  as  with  results 
obtained  in  a previous  study  dealing  with  the  oscillations  of 
the  East  and  West  Basins  of  Long  Beach  Harbor. 

3.  The  coupled-basins  theory  reduces  both  the  computer  time 

and  the  required  computer  storage  in  predicting  the  response 
of  a harbor.  Therefore,  this  increase  in  economy  and 
storage  may  allow  this  type  of  approach  to  be  used  in  pre- 
dicting the  response  of  a harbor  to  shorter  wavelengths  than 
would  be  possible  using  the  single  basin  theory. 

4.  When  harbors  of  simple  geometry  are  connected,  the  response 
of  the  resultant  harbor  system  in  some  ways  appears  to  be  a 
superposition  of  the  response  of  the  individual  harbors.  This 
indicates  that  care  must  be  taken  in  adding  internal  basins  to 
real  harbors  so  that  additional  response  problems  do  not 
occur  after  the  addition  is  completed. 
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LIST  OF  NOTATIONS 


A Wave  amplitude 

A.  Incident  wave  amplitude 

a Characteristic  dimension  of  a harbor,  the  radius  of  a 

circular  harbor 

aA(a_J  Radius  of  the  circular  harbor  referred  in  Basin  A (or 

A a Basin  B) 

B Surface  area 

b Width  of  a rectangular  entrance  channel 

A Dxl  vector  representing  the  values  of  8f21  /8n  at  the 
harbor  entrance,  and  at  the  common  boundaries  between 
Region  II- 1 and  II- 2 as  well  as  between  Region  II- 1 and 
II-3. 

D A number  equal  to  p + dT  + d2 

dj  (d2)  The  total  number  of  segments  at  the  common  boundary 

between  Regions  II- 1 and  II-2  (II- 3) 

Fq  Amplitude  of  an  amplified  force 

Zi  (Za*  Z3)  A Njxl  (N2xl,  N3x1)  vector  representing  the  value  of 

wave  function  at  the  boundary  of  Region  II- 1 (Region  II-2, 
or  Region  II- 3 ) 

f Wave  function  which  describes  the  variation  of  the  velocity 

potential  in  the  x and  y directions 

ft  Wave  function  in  Region  I (in  the  open-sea) 

*ai  (*33  or  *33)  Wave  function  in  Region  H- 1 (Region  II-2  or  Region  II-3) 

f3  Radiated  wave  function 


LIST  OF  NOTATIONS  (Continued) 


Incident  wave  function 


Reflected  wave  function 


An  NjxNi  matrix  defined  in  Eqs.  II,  12e  and  I2f 
(in  Eq.  18  G is  an  NsxNs  matrix  while  in  Eq.  22  G 
is  an  N3xN3  matrix) 

An  NxxNa  matrix  defined  in  Eqs.  11,  12c,  and  12d 
(in  Eq.  18  Gn  is  an  N2xNs  matrix  while  in  Eq.  22  Gn 
is  an  N3xN3  matrix) 


A DxD  matrix  defined  in  Eq.  32 


Zero  and  first  orders  of  the  Hankel  function  of  first  kind 


Water  depth 


Identity  matrix 


Reflection  coefficient 


Wave  number 


Wave  length,  or  representing  a length  scale 
Length  of  a rectangular  entrance  channel 
Logarithm  to  the  Naperian  base  (e  = 2.7128) 


The  ratio  of  the  maximum  excursion  of  the  oscillating 
mass  to  the  static  movement  of  the  mass. 
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LIST  OF  NOTATIONS  (Continued) 


M^  A DxD  matrix  defined  in  Eq.  32 

An  NxxD  matrix  defined  in  Eq.  15 

M2  An  N2xdx  matrix  defined  in  Eq.  20 

An  N3xd2  matrix  defined  in  Eq.  24 

Nj  (N2  or  N^)  Total  number  of  segments  into  which  the  boundary  of 
Region  II- 1 (Region  II-2  or  Region  II-3)  is  divided 

n Outward  normal  to  the  boundary  of  the  region 

Pj  (P2  or  P 3)  A vector  representing  the  value  of  the  normal  derivative 
of  the  wave  function  at  the  boundary  of  Region  II- 1 
(Region  II-2  or  Region  II-3) 

p Total  number  of  segments  into  which  the  harbor 

entrance  is  divided 

R Amplification  factor 

DR  Reynold's  number 

r Distance  between  points  or  radial  position  in  a polar 

coordinates 


Tangent  to  the  boundary  of  the  domain  in  a counter- 
clockwise direction 

Length  of  the  boundary  segments 

Wave  period,  or  representing  a time  scale 
Time 

Is  a NxxD  matrix  defined  in  Eqs.  13  and  14 


Is  a Naxdx  matrix  defined  in  Eq.  19 
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LIST  OF  NOTATIONS  (Continued) 


Is  a N3xd2  matrix  defined  in  Eq.  23 
Velocity  vector  with  components  u,  v,  w 

Coordinate  axis  in  horizontal  direction  parallel  to  the 
coastline 

Position  vector  for  the  point  (x,y) 

Coordinate  axis  in  horizontal  direction  perpendicular  to 
the  coastline 


Function  which  describes  the  variation  of  the  velocity 
potential  in  depthwise  direction  z 

Coordinate  axis  in  vertical  direction 


Quantities  at  the  j segment  of  the  boundary 

The  ratio  of  a quantity  in  the  smaller  model  to  the 
corresponding  quantity  in  a larger  model  or  in  the 
prototype 


Euler's  constant  (y  = 0.577216  ...  ) 

Damping  factor  (see  Eq.  37) 

Displacement  of  water  surface  elevation  from  the  mean 
water  level 


qA(9B) 


Angular  position 

Inclined  central  angle  of  the  opening  of  circular  harbor 
referred  in  Basin  A (Basin  B) 


Exponential  modulus  of  decay  defined  in  Eq.  43 


Circular  wave  frequency  (2tt/T) 


LIST  OF  NOTATIONS  (Continued) 


Velocity  potential 


Forcing  frequency 


Natural  frequency 


Absolute  value 
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APPENDIX 

COMPUTER  PROGRAM 

The  computer  program  for  calculating  the  response  of  an 
arbitrary  shape  harbor  to  the  periodic  incident  waves  which  propagate 
normal  to  the  coastline  is  presented  in  this  Appendix.  In  order  to 
illustrate  the  computer  program,  its  application  to  the  East  and  West 
Basins  of  the  Long  Beach  Harbor  will  be  presented.  Subroutines  used 
in  the  main  program  are  also  listed  for  reference. 

The  main  computer  program  as  well  as  the  subroutines  used  in 
the  main  probram  will  be  listed  first  (pp.  126  to  132).  In  this  program, 
the  region  of  consideration  is  divided  into  three  regions:  an  open-sea 
region.  Region  II- 1 (referred  to  herein  as  Basin  A),  and  Region  H-2 
(referred  to  herein  as  Basin  B).  If  a particular  harbor  geometry  requires 
that  it  be  divided  into  more  regions  then  the  program  must  be  modified  slightly. 
The  input  data  needed  for  using  this  program  are: 

(i)  the  number  of  boundary  segments  of  Basin  A (NA)  and  that 
of  Basin  B (NB),  the  number  of  segments  at  the  entrance 
(NP),  the  number  of  segments  at  the  common  boundary 
between  Basin  A and  Basin  B (ND),  and  the  following  quantities: 
NC  = NP  + ND  The  total  of  the  entrance  segments  plus 
the  total  of  the  segments  at  the  common  boundary  between 
Basin  A and  Basin  B. 


I 


NP1  = NP  + 1 


NM  = The  segment  number  in  Basin  A prior  to  the 
first  segment  of  the  common  boundary  between 
Basin  A and  Basin  B. 

NM1  = NM  + 1 

MA  = The  total  number  of  the  interior  points  in  Basin 
A to  be  calculated. 

NM2  = NM  + ND 

(ii)  the  coordinates  of  the  beginning  and  the  end  of  each  boundary 
segment  in  Basin  B and  Basin  A.  (This  program  is  written 
to  process  first  the  boundary  segments  of  Basin  B and  then 
those  of  Basin  A's;  thus,  it  must  be  supplied  in  this  order.  ) 

(iii)  the  value  of  the  characteristic  dimension  (A),  the  width  of 
the  harbor  opening  (HAOP),  and  the  water  depth  (DEPTH). 
"HAOP"  is  for  identification  only  because  it  does  not  enter 
the  computation. 

(iv)  the  total  number  of  interior  points  in  Basin  A plus  those  of 
Basin  B (M)  and  the  coordinates  of  these  interior  points 
(PX(I),  PY(I)  ),  and 

(v)  the  incident  wave  number  (K) 

These  input  data  for  the  Long  Beach  Harbor  model  for  one 
particular  wave  number  (k  = 2.  35  ft”*)  are  listed  on  p.  133.  There 
are  50  boundary  segments  in  Basin  A and  34  boundary  segments  in 
Basin  B.  The  total  number  of  interior  points  in  Basin  A is  63,  and  the 
total  number  of  interior  points  in  Basin  B is  25.  The  coordinates  of 
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these  boundary  points  and  interior  points  are  arranged  according  to  the 
coordinates  used  in  Fig.  7.  The  number  associated  with  each  interior 
point  listed  herein  is  essentially  the  same  as  used  by  Lee  (1969);  this 
allows  one  to  check  the  results  more  rapidly  in  comparing  both 
methods.  Thus,  point  A of  Fig.  7 corresponds  to  MESH  (26),  point  B 
corresponds  to  MESH  (88),  point  C corresponds  to  MESH  (81),  and 
point  D corresponds  to  MESH  (68). 

The  output  data  for  the  Long  Beach  Harbor  model  are  presented  on 
pp.  134  to  135.  They  contain  the  complex  value  and  the  absolute  value 
of  the  normal  derivative  of  the  wave  function  at  the  harbor  entrance  and 
at  the  common  boundary.  The  complex  value  and  the  absolute  value  of 
the  wave  function  at  the  boundary  of  both  Basin  A and  Basin  B are  also 
shown.  Finally,  the  complex  value  and  the  absolute  value  of  the  wave 
function  (F2)  for  the  88  interior  points  are  printed  on  p.  135;  they  are 
arranged  so  that  those  points  in  Basin  A are  printed  first  followed  by 
those  of  Basin  B.  These  output  results  can  be  checked  with  those 
presented  in  pp.  262  and  263  of  Lee  (1969)  in  which  the  single  basin 
computer  program  was  used  to  calculate  the  response  of  the  same 
harbor  model  and  at  the  same  wave  number  To  calculate  the  response 
at  other  wave  numbers,  data  cards  for  these  wave  numbers  can  be  added 
to  the  last  data  card  presented  here.  It  is  also  noted  that  this  program 
is  written  in  FORTRAN  IV  compatible  with  the  IBM  360/75  digital 
computer. 
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Some  of  the  symbols  used  in  the  computer  program  and  not 


previously  defined  in  this  appendix  are  as  follows: 


M = Total  no.  of  interior  points  to  be  calculated 

(Basin  A plus  Basin  B) 

NSEG(I)  = Number  which  defines  the  boundary  segments  of 

Basin  A or  Basin  B (also  used  to  define  the  interior 
points  to  be  calculated) 

fch 

PX(I)  = The  x-coordinate  at  the  beginning  of  the  1 segment 

of  the  boundary  of  Basin  A (or  Basin  B)  (also  used 
as  the  x-coordinate  of  the  interior  points) 

til 

PY(I)  = The  y-coordinate  at  the  beginning  of  the  i segment 

of  the  boundary  of  Basin  A (or  Basin  B)  (also  used 
as  the  y-coordinate  of  the  interior  points) 

MESH(I)  = Number  which  defines  a particular  interior  point 


EKA 

R(I,  J) 

DX(I)  (or 
DY(I)  ) 


DS(I) 


PERT 


Wave  number  parameter  (ka) 

Distance  between  field  and  source  points 

x (or  y)  projection  of  the  length  of  the  i n boundary 
segment  of  Basin  A or  Basin  B 

Length  of  the  i^  boundary  segment  of  Basin  A or  Basin  B 
Wave  period 


DFDN  (I,  J)  = An  Nj  x D matrix  (for  Basin  A)  equivalent  to  the  matrix 

defined  in  Eqs.  13  and  14;  for  Basin  B it  is  an  N2xdj  matrix 


AN  (I.  J)  = 


DRDN 


XSS  (YSS)  = 


TEMP 


equivalent  to  matrix  defined  in  Eq.  19. 

A matrix  defined  as  (yG  - l\  it  represents  an  N,xN, 

\2  n r 11 

matrix  for  Basin  A,  and  is  an  N2  x N2  matrix  for  Basin  B. 
9r/9n 


9*x/9s3  (93y/9s3 ) 

As  / 9x  93y  9ax  9y  '\ 
it  \9s  9sa  " 9s3  9s  / 


| I 
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G(I.  J) 


An  N,xNj  matrix  for  Basin  A (Eqs.  12e  and  12f)  or  an 
^xN^  matrix  for  Basin  B (Eq.  18) 


Q(I.  J) 


An  N^xD  matrix  for  Basin  A equal  to  the  matrix  ^-GU^ 
defined  in  Eq.  15;  for  Basin  B it  is  an  ^xd^  matrix 


CSLECD  = 


A subroutine  for  solving  complex  systems  of  linear 
equations 


Q(I,  J) 


(After  the  statement  CALL  CSLECD)  It  represents  the 


matrix  Mj  defined  in  Eq.  15  for  Basin  A;  for  Basin  B it 
represents  M2  defined  in  Eq.  20. 


A DxD  matrix  defined  in  the  process  of  matching  (see 
Eq.  33) 


DFDC  (I.  1)  = 


A Dxl  complex  numbered  vector  represents  the  normal 
derivative  of  the  wave  function  at  the  entrance  and  that 
at  the  common  boundary 


ADFDC  (I,  1 ) = Absolute  value  of  DFDC  (I,  1) 


Q(I.  1) 


Complex  number  representing  the  value  of  the  wave 
function  at  the  boundary  segments  of  Basin  A 


QB  (I,  1)  = 


Complex  number  representing  the  value  of  the  wave 
function  at  the  boundary  segments  of  Basin  B 


ABBF  (I,  1)  = Absolute  value  of  Q (I,  1)  (or  QB  (I,  1)) 


FI  (I) 


The  complex  value  of  f for  the  interior  points  (either 
in  Basin  A or  in  Basin  B) 


F2  (I) 


Absolute  value  of  FI  (I)  with  the  sign  equal  to  that  of  the 
real  part  of  FI  (I) 


FR  (I) 


The  ratio  of  F2  (I)/F2  MAX 


MAXMIN 


Subroutine  to  find  maximum  and  minimum  elements  of  an 
array 


INTEGER  P 
REAL  K 

COMPLEX  AN( 50, 50 ) «Q( 50, 10 ) »QB( 34,10)«H(50.10) ,0F0C(50,3 ) ,G(50.50 ) 
COMPLEX  Fl( 100) .OET.C.D.OGON.F 

DIMENSION  AOFDCI 50.1) ,DCDN (50.10) . ABBF ( 50. 1 ) *F2 ( 100  I *FR( 100), 

1 PX( 100) ,PV(100) ,NSEG( ,00) ,R (50.50 > ,RA ( 50.50 ) ,RB( 34 ,34 ) • 

2 X0(1),X(51).Y0(1).Y(51).DS0(1),DS(51).DX(50).DY(50). 

3 BX0(1).BX(35) ,BYO< l) ,BY(35) ,DSB0(1 ).DSB(35),DXR(34),DY8<34). 

4 AXO( 1 >,AX(51 ).AY0(1).AY(51),DSA0( 1) .DS4(5l) ,0X A ( 50 ) ,OYA ( 50 ) 

DATA  PI/3.1415926/ 

PI3*3.0*PI 

C«CMPLX(0.,-l. 0/2.0) 

D*CMPLX (0. ,-0.25 ) 

THOOP I *2 .0/P I 
C READ  INPUT  DATA 

READ  (5,1)  NA,NB,NP,ND,NC.NPl,NMfNMl,MA,NM2 

1 FORMAT  (1415) 

N*NB 

KT I ME* 1 

11  READ  (5,2)  (NS6G( I ) ,PX( I ) ,PY( I ) ,1=1 ,N> 

2 FORMAT  I15.2F13.0) 

NhR  = N - 1 

px(N41)*px( n 
PY ( N+ 1 ) = PY ( 1 ) 

PX ( N>2 ) *PX ( 2 ) 

PY ( N^2 ) =PY ( 2 ) 

C CALCULATE  MIDPOINT  OF  EACH  SEGMENT 

00  5 1*1, N 

X(  I )*0.5*(PX(  I )*PX(  I + l ) ) 

Y(  I )*0.5*(PV(  I ) ♦PY  ( 1+1)) 

DX( I )*PX( I *1 )-PX( I ) 

DY ( 1 )*PY( I ♦ 1 ) — PY ( I ) 

0S( I )*SORT(OX( I ) ♦*2+DY ( I ) **2 ) 

5  CONTINUE 
X0( 1)=X(N) 

X ( N+l ) *X ( 1) 

Y0( 1)*Y(N) 

Y ( N+  l ) * Y ( 1) 

0S0(1 )*0S(N> 

D$( N+ 1 ) = DS ( 1 ) 

R ( N.N)  * 0 
DO  15  1*1 »NMR 
Il»IM 
R ( I • I ) *0 
DO  25  J* I l , N 

R ( I ,J)*SORT( (X( I )-X(J) ) ♦*2* ( Y ( I ) -Y ( J ) ) **2 ) 

R ( J , I * *R ( I ,J) 

25  CONTINUE 
15  CONTINUE 

IF  (XTIME.E0.2)  GO  TO  100 
C STORE  HARBOR  GEOMETRY  OF  THE  SECONARY  BASIN 

BXINB-U  )*X(NBd  ) 

BY(NB+1)*Y(NB*1) 

DSB ( NB* 1 ) *0S ( NB* 1 ) 

DSBO ( 1 ) *DSO( 1 ) 

8Y0(1)*Y0(1) 

BXO ( l)*XO( 1) 

DO  415  1*1. NB 
OXR( I ) *DX ( I ) 

DYB ( I )*DY( I ) 

DSB ( I ) *DS( I ) 

BX( I )*X(  I ) 

Bv ( I ) *y( I ) 

DO  415  J* 1 , NB 
R B( I • J ) *R ( I.J) 

415  CONTINUE 
N*NA 

XT I ME  * 2 
GO  TO  11 
100  CONTINUE 

STORE  GEOMETRICAL  INFORMATIONS  OF  FIRST  BASIN 
AX(NA*1 )«X(NA*1 ) 

AY(NA*1)*Y(NA*1) 

DSA ( NA* 1 ) *0S  t NA* 1 ) 

AXOll)-XO(l) 

AY0I1)*Y0U) 

OSAO(l)-OSOd) 

DO  425  I«1 ,NA 
AX(|)-X(|) 

AYII)-Y(I) 

OXA ( I ) *0X(  I ) 

DVAI I )*DV( | ) 

OSA(I)-OSII) 

00  425  J-l.NA 
RAII, J)-R( I.J) 

425  CONTINUE 


C 


upryi  '7*~' ■ 
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REAO  (5*4)  A 
4 FORMAT  (4F10.0) 

REAO  (5*4)  HAOP, DEPTH 

REAO  COORDINATES  OF  INTERIOR  POINT  INTO  PX  AND  PY 
READ  (5,1)  M 

READ  (5.2)  (NSEG( I),PX( I),PY(I),I*l,M) 

16  READ  (5.17)  K 

17  FORMAT  (FIO.O) 

EKA»A*K 

PERT*(2.0*PI >/(SORT(32.2*K*TANH(K*DEPTH) ) ) 

CALCULATE  UNIT  MATRIX  DFDN  FOR  SECONARY  BASIN 

N»NB 

NN*ND 

XO( 1 )«BXO( 1 ) 

YO ( 1 ) *BYO( 1) 

DSO( l)-DSBO( l) 

X ( NB*1 )*BX(NB+1 ) 

VINBMl-BYlNB+l) 

DSINR+l  )*DSB(  NB*  I ) 

DO  435  1*1, NB 
X ( I ) *BX  ( I ) 

Y ( I ) *BY I I ) 

DXU >=DXB( I ) 

OY( I ) *DYB ( I ) 

OS< I )=DSB( I ) 

00  435  J=1,NB 
R ( I . J I *R8( I » J ) 

435  CONTINUE 
110  CONTINUE 

DO  175  I*l.NB 
DO  175  J=1,N0 

IF  ( (I+J)-(ND+1>  ) 185,190.185 
190  DFDN( I ,J)*-1.0 
GO  TO  175 
185  DFDN( I ,J)=0.0 
175  CONTINUE 
GO  TO  60 
160  CONTINUE 

CALCULATE  UNIT  MATRIX  DFDN  FOR  MAIN  BASIN 
DO  115  1=1, NA 
DO  115  J=ltNC 
OF ONI  I , J ) =0.0 
115  CONTINUE 

DO  12*  1=1, NP 
OFONU, II*  1.0 
125  CONTINUE 

DO  105  J=NP 1 , NC 
1 =NM-NP+J 
DFDNI I , J ) = 1.0 
105  CONTINUE 
60  CONTINUE 

CALCULATE  ELEMENTS  OF  THE  MATRIX  AN=(C*GN-I) 

DO  35  1*1  ,N 
DO  45  J=l ,N 
IF  ( J .EO.  I ) GO  TO  10 
ARG*K*R< I « J ) 

DRDN* ( ( Y ( I I -Y ( J ) ) *DX ( J ) - ( X ( I ) -X ( J ) ) *DY I J ) >/R( I ,J ) 

AN( I « J ) =-C*K*CMPL  X ( BE  S J 1 ( ARG ) , BE  SY 1 ( ARG ) )*DRON 
GO  TO  45 
10  CONTINUE 
YSS= 

*6.0*( (VI  I ♦ 1 ) — Y ( I) I / ( D S I I + 1)+DS< I ) ) - ( Y ( I ) — Y ( I — 1 ) > ✓ ( 0 S ( I ) ♦0  S ( I-l)  ))/ 

/ (DS(  i-n+osm+osi  IM)  ) 

XSS* 

*6.o*(  ( xi  iM)-x(i))/iDSiiM)*Dsm)-<xiii-xii-in/iDsin+Dsii-im/ 
/ (DS( I-l )*DS( I )*DSl 1*1 ) I 
TEMP* | DX I I )*Y$S-XSS*DV< I ) )/Pl 
AN( I, I >=CMPLX(0.0,TEMP)*C-1.0 
45  CONTINUE 
35  CONTINUE 

CALCULATE  THE  RIGHT  HANO  SI06  VECTOR  U 
DO  135  P*1,NN 
DO  55  1*1, N 
Q<  1 ,P)*CMPLX(0.,0. ) 

DO  65  J* 1 • N 
IF  (P.NE.l)  GO  TO  30 
IF  I J .EO.  I ) GO  TO  20 
ARG*K *R I I ,J> 

G(  I , J )*CMPLX ( BES J0( ARG) ,BE SVO ( ARG ) ) *DS  ( J ) 

GO  TO  30 

20  G( 1,1 )*CMPLX( 1.,TmOOPI*(ALOG(K*DS( J)*.25)-0. 42279) ) *DS ( I ) 

30  0(  I ,P|«0(  !,P)+GU  ,JI *OFON<  J,P) 

65  CONTINUE 

0(  I.P)«C*0( l,P) 

55  CONTINUE 
135  CONTINUE 

CALL  CSLECDI AN,N,0,NN,DET,IER) 

IF  (N.EO.NA)  GO  TO  70 


STORE  INFORMATION  FOR  SECONOARV  BASIN 
DO  255  I * 1 » NB 
00  255  J * 1 • NO 
OBI  I * J) *01 1 » J) 

255  CONTINUE 
N*NA 
NN*NC 

xom-Axoii) 
vom«AYom 
OSOI 1 ) *DSA0l I ) 

X(NA+1)*AXINA+1) 

Y I NA+1 )*AY(NA+1) 

DSINA+1 )=DSA(NA+1) 

00  275  1*1 , NA 
X | I ) *AX ( I ) 

Yl  n*AYI  I ) 
oxm*oxAi  i ) 

OYI  n*0YA(  I ) 

OSI I )=DSA| 1 ) 

DO  275  J*l,NA 
R« I * J ) *RA(  I ♦ J I 
275  CONTINUE 
GO  TO  160 
70  CONTINUE 

00  155  I * 1 * NA 
00  155  J*l,3 

OFDCI I, J)=CMPLX(0.0,0.0) 

155  CONTINUE 

CALCULATE  WAVE  FUNCTION  OF  EXTERIOR  PROBLEM  ANO  MATCHING 
00  205  1 = 1 * NC 
00  205  J=1 ,NC 
HI  I » J ) =CMPLX 1 0.0* 0.0 ) 

205  CONTINUE 

00  215  1*1 » NP 

DO  215  J= 1 * NP 
IF  I I.EQ.J)  GO  TO  210 

HI  I , J ) *CMPLX( BE  S JO  I K*R I I , J) > ,HESYO(K*R( I , J ) ) >*DS( J )*C 
GO  TO  215 

210  HI  I « I ) *CMPLX I 1 ,0«T WOOP I * I A LOG  I K*DS( I ) * . 25 ) -0.52279 ) )*DS( I ) *C 
215  CONTINUE 

00  315  1 =NP l » NC 
DO  315  J=NP 1 ,NC 

1 I =NC- !♦ 1 
JJ= J-NP 

HI  I * J ) *0B( 1 1 * J J ) 

315  CONTINUE 

DO  325  1*1 »NP 

00  325  J * 1 * NC 

HI  I t J ) *0  I I , J ) -*H( I *J) 

325  CONTINUE 

DO  335  I =NP1 * NC 
DO  335  J* 1 » NC 

1 I *NM+ I-NP 

HI  I , J ) *01  1 1 , J ) "Hi  I , J ) 

335  CONTINUE 

DO  225  1*1* NP 
OFDCI I *1 )=CMPLX< 1.0,0.01 
225  CONTINUE 

CALL  CSLECD(H,NC,DF0C,1,0ET,!ER) 

WRITE  16,6) 

6 FORMAT  I 1H1) 

WRITE  16,19)  K 

19  FORMAT  I2X3HK  *F  10.5 , 2X, • I 1/FT  ) • ) 

WRITE  16,38) 

38  FORMAT  I / // , 2X , • COMPLE X VALUE  OF  DFOC  AT  THE  ENTRANCE  ANO  AT  THE  C 
♦ OMMON  BOUNOARY  I I /FT)*,/) 

WRITE  16,8)  (OFDCI 1,11, 1*1, NC) 

8 FORMAT  I IX, 6F13 .5) 

00  305  I * 1 , NC 

305  AOFDCI 1,1 )*CA8S(0F0CI 1,1)) 

WRITE  (6,68) 

68  F0RMATI///,2X,» ABSOLUTE  VALUE  OF  DFOC  AT  THE  ENTRANCE  11/FT)*,/) 
WRITE  (6,8)  I AO FOCI I,1),I*I,NP) 

WRITE  16,69) 

69  FORMAT  (///, 2X ,* ABSOLUTE  VALUE  OF  DFOC  AT  THE  COMMON  BOUNDARY  1 1 /F 

•TIM 

WRITE  (6,6)  (AOFDCI 1,1 ),I*NP1,NC) 

CALCULATE  BOUNOARY  WAVE  FUNCTION  F OF  PRIMARY  BASIN 
00  235  J-I,NA 
DO  235  J*  1 • NC 

235  OFDCI I ,2 )*DF0CI I ,2)+OI I , J)*OFDCI J* 1 ) 

DO  255  1*1, NA 
011,1 )*0F0CI 1,2) 

ABBFI t,l)-CABSIO( 1,1) ) 

255  CONTINUE 


annca 
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WRITE  (6*48) 

48  FORMAT  (///* 2X *• COMPLEX  VALUE  OF  THE  BOUNDARY  F FUNCTION  0(1*1)*) 

WRITE  (6*47) 

47  FORMAT  (2X, '(BASIN  A)*,/) 

WRITE  (6*8)  (0( 1*1) ,1*1. N) 

WRITE  (6*148) 

140  FORMAT  (///*2X*' ABSOLUTE  VALUE  OF  THE  BOUNDARY  F FUNCTION') 

WRITE  (6*47) 

WRITE  (6*8)(ABBF( 1*1),1*1,N) 

C CALCULATE  BOUNDARY  WAVE  FUNCTION  F OF  SECONDARY  BASIN 

DO  385  1*1, NB 
UFUCI I ,2)*CMPLX(0.0.0.0) 

OF DC ( I*3)»CMPLX(0.0*0.0) 

DO  375  J * 1 , NO 
JJ*NP>J 

375  OFDC(  1*2)  *DFOC(  1*2)  ♦QH(  I , J ) *OFI)C ( J J , 1 ) 

0B< 1 • 1 )*DFDC( I *2) 

ABBF ( I , 1 ) *CAHS ( Ortf 1,1)) 

385  CONTINUE 

WRITE  (6*48) 

WRITE  (6*49) 

49  FORMAT  I2X, '(BASIN  8)',/) 

WR I T F ( 6*  8 ) ( OH ( I • 1 ) • 1 * 1 * N8 ) 

WRITE! 6* 148) 

WRITE  (6.49) 

WR I TF ( 6*  8 ) ( AHBF(  I , l > , 1*1 ,NH) 

C NUMBERING  BOUNDARY  F DERIVATIVES  FOR  SECONDARY  BASIN 
DO  355  1*1, NO 
L=NC- I ♦ l 

DFOCI  1*3)  =DF()C  ( L , 1 ) 

356  CONTINUE 

C RENUMBERING  BOUNOARY  DERIVATIVES  OF  PRIMARY  BASIN 

DO  365  I *NM l * NM 2 
J* I -NMl+NP 1 
OFOC< l « 1 1 *0F0C( J , 1 ) 

DFDC( J, 1 )*CMPLX(0.0,0.0) 

365  CONTINUE  j 

C CALCULATE  WAVE  FUNCTION  F FOR  INTERIOR  POINTS 

DO  75  J=1,M 
F*CMPLX(0. 0,0.0) 

IF  (J.GT.MA)  GO  TO  395 
DO  85  I * 1 « NA 

Rl* SORT ( ( AX  I I ) -PX ( J ) ) **2+  < A Y( 1 ) -RY( J ) ) **2 ) 

RK«K*R 1 

DGON=K*CMPL X ( HE  S J 1 ( RK ) , BE S Y l ( RK ) ) * ( ( PX ( J ) -AX ( I ) )*DYA( 1 ) 

* -IPY( J)-AY(  I ) )*DXA(  I ) )/Rl 

F*F*0( I , 1 )*DGON-DFDC( I . I ) *CMPlX ( BE SJO ( RK ) * BE SYO ( RK  ) )*OSA( I ) 

85  CONTINUE 
GO  TO  200 
395  DO  495  1*1, NB 

Rl * SORT ( (BX(  I )-PX(J) ) **2+ ( BY ( I ) -PY ( J ) ) **2  ) 

RK*K*R1 

DGDN*K*CMPLX ( BE  SJ 1 ( RK ) *BESY1 ( RK ) ) *(  (PX( J)-BX( I ) )*OYB( I ) 

♦ -( PY ( J )-8Y (I))*0XB(I))/R1 

F*F  + OB ( 1*1) •OGON+DFUC ( I * 3 ) *CMPL X ( 0E S JO ( RK I * BE SYO < RK  I )*OSfl ( I ) 

495  CONTINUE 
200  F*D*F 
Fil J)*F 

F2( J ) *S I GN( CABS ( F ) ,REAL(F) I 
75  CONTINUE 

CALL  MAXMIN( F2*M*F2MX*F2MN ) 

F2MAX*AMAX1( ABS( F2M X ) * ABS ( F 2MN ) ) 

DO  605  J* 1 *M 
605  FR(J)*F2( J)/F2MAX 
WRITE  (6*6) 

WRITE  (6*119)  HAOP.OEPTH 

119  FORMAT  (2X, 'HARBOR  OPENING  (FT.)«',F  7.3.5X  , 'DEPTH  (FT.)*'«F7.3) 

WRITE  (6*19)  K 
WRITE  (6*26)  EKA 
26  FORMAT  ( 2X3HKA*F 10.5 ) 

WRITE  (6*129)  PERT 

129  FORMAT  (2X*'PERI00  T *• ,F 10.5 , 2X *•( SEC. >•  ) 

WRITE  (6*199)  F2MAX 
199  FORMAT  ( // , 2X • • F2MAX- • , F 10.5 ) 

WRITE  (6*99) 

99  FORMAT  ( // * 2X , • ME SH* * BX.'PX',  9X*'PY' • 14X  * 'F CMP LX  *«15X«  'F2»* 

• l OX  * • FRA • • , / ) 

WRITE  (6*9)  (NSEG( J),PX( J ) ,PY ( J ) ,F 1 ( J ) ,F2 ( J)*FR( J)*J-1,M) 

9 FORMAT  (U*I5*F11.3*FU.3*2F13.5*F13.5*F11.3) 

GO  TO  16 
ENO 
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SUBROUTINES 


(1) 


SUBROUTINE  CSLECDI  A * M,  B • N • OET,  IUI 
C SOLUTION  OF  COMPLEX  SYSTEM  OF  L IN.EOUAT .W I TH  N RIGHT  HAND  VECTORS 

C ANO/OR  COMPUTATIUN  OF  COMPLEX  DETERMINANT 

COMPLEX  A * H,  AT,  FAC.  OET 
DIMENSION  A( 50  »M) , B(50«N) 

ILL*  0 

CALL  OVERFLI 10) 

SIGN*  ♦ ! 

IMA*  M- 1 
DO  35  1*1, IMA 

AMAX*  REAL  < A ( 1,1  ) ) « RF  AL I A < | « I » ) ♦ A I M A G ( A < I , I ) ) * A I M AG ( A ( I , 1 ) ) 

JMAX*  1 

II*  U1 

00  20  J*  1 1 ,M 

ARE*  REAL I AI  J , I \ > 

AlMs  A I MAG( A ( J , I ) ) 

A J I = ARE*AR8  ♦ AIM*AIM 
I F ( AMAX-A J 1 ) 18.20,20 

18  AMAX*  A J I 
JMAX*  J 

20  CONTINUE 
IF(AMAX)  21.90.21 

21  IF! I — JMAX ) 23,25.23 

23  SIGN*  -SIGN 
00  2*  K* I ,M 
AT*  All .K) 

A | I ,K ) * AlJMAX.K) 

24  AIJMAX.K)*  AT 
IFIN.LE.O)  GO  TO  25 
DO  241  K * 1 « N 

AT*  R! 1 ,K) 

B(I.K)*  BIJMAX.K) 

241  BIJMAX.K)  * AT 

25  DO  35  J*  1 1 .M 

FAC*  A(  J.  H/A(  Il-I  . I l-l  I 
DO  30  K*U,M 

30  A ( J • K ) * A I J, K ) - FAC*A I I , K ) 

IFIN.LE.O)  GO  TO  35 
DO  32  K*1,N 

32  BIJ.K)  = BIJ.K)  -FAC*BU.K) 

35  CONTINUE 

C TRIANGULAR  MATRIX  RFADY 

IFIN.LE.O)  GO  TO  70 

,F  ICABSI AIM,*) ) .EO.  0.)  GO  TO  90 

DO  40  K*1.N 

40  8(M,K)  * 8 I M , K ) / AIM.M) 

DO  60  1*1,  IMA 

J*  M-I 
K 1 = JM 
DO  50  K*K  1 ,M 
DO  50  L*l.N 

50  BIJ.D*  8 ( J , L ) -A!  J,K)*BIK,L) 

DO  60  L*l.N 

60  BIJ.D  * BIJ.D  / AIJ.J) 

70  OET*  A 1 1 « 1 ) 

DO  74  1*2, M 
74  DET*  DET*A I 1,1) 

OET*  DET*  SIGN 
CALL  OVERFLI 10) 

IFIIO.EO.D  GO  TO  91 
RETURN 

90  OET*  10,, 0. ) 

91  WRITE! 6,92 ) 

92  FORMAT  1 46H0DET  A * 0 OR  OVERFLOW  IN  SUBROUTINE  CSLECD  ) 

ILL*  -1 

RETURN 

END 


FUNCTION  BESJOIX) 

BESSEL  FUNCTION  JO(X) 

IF  (X  .GT.  3.)  GO  TO  10 
IF  (X  .LT.-3.)  GO  TO  20 
X03-X/3. 

T«X03*X03 

BE S JO* l • ♦T*C-2.24999974TI*| 1 .26562084T ♦( -.3 163866  ^T*! .0444479 
1 ♦T*f-.  0039444  ♦T*.0002100l ) ) I I 

RETURN 
10  T-3./X 

F0«.79788456^T*(-.00000077^T*( - .0055 2740^T*(- .000095 12 
l 4T*I .OOl37237^T*l -• 000 72805 ♦T*. 000 144 76 ) ) ) ) ) 
THETA0«X-.78539816^T*(-.041663974T* I -.00003954 
1 ♦T*! »00262573^T*I -.00054125^T*I -• 000293 3 3^T 4.000 1355 8 ) I I » 1 

8ESJO*FO*COSf f HE  T AO  > / SOR  T I X ) 

RETURN 

20  WRITE  (6.6) 

6 FORMAT  I44H  ARGUMENT  LESS  THAN  -3  * NO  JO(X)  CALCULATED) 

RETURN 

END 


FUNCTION  BESJ1IX) 

BESSEL  FUNCTION  J 1 1 X » 

IF  IX  .GT.  3.)  GO  TO  10 
IF  CX  .LT.-3. ) GO  TO  20 
X03*X/3. 

T*X03*X03 

BE  SJ 1 = • 5^T* I - • 5 62499 85 ♦T  * ( .2 l 0935 73 ♦T  * I - .039542B9^T ♦ I .00443  319 
1 ^T* I -.0003 l 761 ♦.00001 109) ) ) ) ) 

BESJl*BESJl*X 
RETURN 
10  TO./X 

Fl*.79788456^T*f .000001 56+T ♦ I .01659667^T4| .00017105 
1 ♦T*  (-.  002495 1UT*1  .001 13653-T».  00020033)  ) ) )) 

THETA1 *X-2.35619449^T  ♦ ! . 124996124T*! • 00005 650^T* I -.00637879 
1 ♦ T*! .000743484T*! .00079824-T* .00029166 ) ))) ) 

BE SJ 1* F 1 *COS ( THE  T A 1 ) / SORT  IX) 

RETURN 

20  WRITE  16.6) 

6 FORMAT  I 44H  ARGUMENT  LESS  THAN  -3  • NO  Jl(X)  CALCULATED) 

RETURN 

END 


n 


FUNCTION  BESYOIX) 

NEUMANN  FUNCTION  YO(X) 

IF  IX  .LE.  0. ) GO  TO  20 
IF  IX  .GT.  3.)  GO  TO  10 
X03*X/3. 

T=X03*X03 

BESY0  = 2./3. 14 159  2 6*  A LOG  I . 5*X  ) ♦BE  S JO  I X ) ♦ . 3674669UT  ♦ I .60559366 

1 ♦T* I - • 743  503  84  + T ♦ I . 25300 1 1 7^T ♦ I - .0426 1 2 14^T ♦ ( .00427916 

2 -T*. 00024846) ))) ) 

RETURN 

10  T*3. /X 

FO* • 7978845 6*T* I - .000000 77 ♦T * ( - .0055  274 O^T ♦ I - .000095 12 
1 ♦ T4 1 • 00 137237^T*I— * 000  72805^T*. 000 14476) ) ) ) ) 

THETA0*X-.78539816^T  *1  - .04 1663 97^T* (- .00003954 
1 ♦T*! . 002  62573>T^|-#  00054 125^T*I-.00029333^T4. 00013558 ) ) ) ) ) 

BESY0*F0*SINlTHETA0) /SORT ( X ) 

RETURN 

20  WRITE  16.6) 

6 FORMAT  (48H  ARGUMENT  LESS  THAN  OR  * 0 • NO  YO(X)  CALCULATED) 
RETURN 
ENO 


FUNCTION  BESV1IX) 

NEUMANN  FUNCTION  Yl(X) 

IF  IX  .LE.  0.)  GO  TO  20 
IF  IX  .GT.  3.)  GO  TO  10 
X03-X/3. 

T*X03*X03 

BESYI*2./3.1415926<X4ALOGI0.54X)*BESJ1IX) 

1 -• 6366198^T ♦! .2212091^T4«2.1682709>T*I-1 .3164827 

2 ♦T*! • 3l2395l*T*l-.0400976^T*. 0027873) ) ) ) ) 

BESV1-BESY1/X 

RETURN 
10  T-3./X 

F1».797884564T*I .000001564T4I .016596674T*! .00017105 
l < -.0024951 1 ♦T • I .001 1 3653-T*. 00020033 ) ) ) ) ) 

THETA1-X-2. 3561944 9^T*t . 1249961 2*T*I ,00005650+T ♦ I - . 00637879 
l ♦T*I.0007434B4T*(.00079824-T4. 00029166) ) ) ) ) 

BESVI«F1*S1NITHETA1) /SORT  I X ) 

RETURN 

20  WRITE  (6*6) 

6 FORMAT  ( 48H  ARGUMENT  LESS  THAN  OR  ■ 0 « NO  Yl(X)  CALCULATED) 
RETURN 
END 


THE  SUM  OUT  INC  •NAXNIN*  F|NO  S 7H€  UMe*  AMO  KMM  MHJMDAA 16  S 
OF  AN  AAA  AY  * CITHER  R€  AL  OR  INTEGER.  THE  CALLING  SEQUENCE  IS 
AS  FOLLOWS* 


WHERE  A - ONE  DIMENSIONAL.  AAAAV*  REAL  OA  INTEGRA* 

N • LENGTH  OF  THE  AAAAV.  N MUST  BE  CREATE A OR  EOMAL  TO  I 
• POSITIVE*  IF  A IS  A FLOATING-POINT  AAAAV. 

■ NEGATIVE*  IF  A IS  AN  INTEGER  AAAAV. 

ANX  * MAX  I HUH  VALUE  OF  THE  ARRAY*  REAL  OR  INTEGER. 

ANN  ■ N IN (HUM  VALUE  OF  THE  ARRAY*  REAL  OR  INTEGER. 


Au  EOli  0 

R1  EOU  I 

A 2 EOU  ? 

A3  EOU  3 

I EOU  A 

A EOU  S 

ANX  EOU  6 

ANN  EOU  7 

BASE  EOU  8 

N EOU  9 

A ANX  EOU  10 

AANN  EOU  1 i 

FNX  EOU  0 

FNN  EOU  2 

FR  EOU  A 

ENTRY  NAXHIN 

NAXNIN  SAVE  U4*l?)**« 

BALA  BASE  *0 

USING  ••BASE 

L A AH  X * 8 ( R 1 ) GET  AOOR.  (IF  AMX. 

L AAMN*I?(A1)  GET  AODR*  UF  ANN. 

L R2*4| All 

L N*0IR2)  GET  N. 

C N* *F * | • IS  N * | 7 

BNF  NNE I NO*  BRANCH. 

L A2  *0 1 R 1 ) YES.  SET 

L A2.0IR2)  AMX  ■ All)* 

ST  R2  * 0 1 AAH  X ) AMN  « AU). 

ST  A 2* 01 AAHN ) AND  RET  URN  • 

RE TN  RETURN  (14*|2)*T 

NNE 1 LTR  N*N 

BNZ  NNZ  BRANCH  |F  N *NE • 0. 

L 1S.»V( IBCOHAI  IF  N * 0*  PRINT 

CNOP  0*4  AN  ERROR  HFSSAGt 

HAL  1 4 , 4 ( IS)  ANO  TERMINATE  THE  JOB 

DC  A(6)  *AL1  IM  *AL3(FM) 

BAL  14.I6I1S) 

CALL  EXIT 

NNZ  L A.OIRl) 

LPR  R3.N 

BCTR  R 3 * 0 

SLL  R3  *2 

LA  1*4 

LA  R2*4 

LTR  N*N 

BP  FLOATPT 

L ANX, 0(A) 

LR  ANN* ANX 

LOUP  1 L R0*0( A* I ) 

CR  ANX  *R0 

BL  SETAHX 

CR  AMN* HO 

bnh  bxle 

LR  ANN*  RO 

B BXLE 

SET ANX  LR  ANX  *R0 

BXLE  BXLE  l*R2*LU0PI 

ST  ANX  OIAANX) 

ST  ANN*  0 1 A AMN ) 

B RETN 

FLOATPT  LE  FNX*0(AI  SET  FNX  * Alii* 

LE  FNN, OCA)  FNN  * AI1)  INITIALLY* 

L00P2  LE  FR*0<A*!>  GET  AIM. 

CER  FNX  * FR  IS  FNX  >«  All)  ? 

BL  SET  FNX  NO.  BRANCH. 

CER  FNNfFR  IS  FNN  <«  All)  7 

BNH  FBXLE  YES.  BRANCH. 

LER  FNN.FR  NO.  SET  FNN  •All) 

B FBXLE  AND  BRANCH. 

SETFNX  LER  FNX.FR  IF  FNX  < All),  SET  FNX  •All) 

FBXLE  BXLE  l*R2*LOOP2  STEP  1 AND  LOOP  BACK. 

STE  FNX, OIAANX)  STORE  THE  FLOATING-POINT 

STE  FNN* 01 AANN ) RESULTS  ANO  RETURN 

B RETN 

FN  OC  CM/"  ERROR  RETURN  FRON  NAXNIN— N ■ 0.*»)» 

ENO 


initialize  ARRAY  POINTER. 

REG. 3 « IABSIN). 

SET  UP  UPPER  UNIT  OF  ARRAY 
BY  SETTING  REG. 3 * (N-l)*4. 
SET  I * 4 TO  START  THE  LOOP. 
REG. 2 * 4,  INCREMENT  UF  I. 

IS  N PLUS  ? 

YES.  A IS  FLOAT ING-PT.  BRANCH 
SET  AMX  « AU)  • 

AMN  * All)  INITIALLY. 

GET  AU). 

IS  AMX  >«  AIM  7 
NO*  BRANCH* 

IS  AMN  <•  All)  7 
YES*  BRANCH* 

NO.  SET  AMN  * AU) 

ANU  BRANCH. 

IF  AMX  < AU  ) , SET  AMX  • AU  I 
STEP  1 ANO  LOOP  BACK. 

STORE  THE  INTEGER  RESULTS 
ANO  RETURN. 


INPUT  DATA 
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88 


1 

-0.210 

-0.106 

2 

0.750 

-0.225 

\ so 

34  2 

4 6 

3 

7 

8 

63 

11 

3 

0.600 

-0.225 

1 

-0.285 

-0.698 

4 

0.450 

-0.225 

2 

-0.330 

-0.64-* 

5 

0.300 

-0.226 

3 

-0.390 

-0.570 

6 

0.150 

-0.225 

4 

-0.450 

-0.496 

7 

0.000 

-0.225 

S 

-0.491 

-0.450 

8 

-0.150 

-0.226 

* 6 

-0.SSS 

-0.4  66 

9 

-0.270 

-0.226 

i 7 

-0.675 

-0.486 

10 

1.035 

-0.405 

8 

-0.809 

-0.509 

1 1 

0.900 

-0.376 

},  9 

-0.930 

-0.530 

12 

0.750 

-0.375 

? 10 

-0.990 

-0.545 

13 

0.600 

-0.375 

1 1 

-1.023 

-0.585 

14 

0.450 

-0.375 

12 

-1.073 

-0.720 

15 

0.300 

-0.3  76 

13 

-1.115 

-0.842 

16 

0.150 

-0.375 

14 

-1.161 

-0.975 

17 

0.000 

-0.375 

IS 

-1.199 

-1.079 

18 

-0.150 

-0.375 

16 

-1.202 

-1.115 

19 

-0.300 

-0.376 

17 

-1.179 

-1.146 

20 

1.200 

-0.525 

18 

-1.079 

-1.184 

21 

1 .060 

-0.525 

-0.954 

-1.229 

22 

0.900 

-0.526 

20 

-0.825 

-1.280 

23 

0.760 

-0.526 

21 

-0.774 

-1.302 

24 

0.600 

-0.525 

22 

-0.719 

-1.310 

25 

0.450 

-0.526 

23 

-0.686 

-1.335 

26 

0.300 

-0.525 

24 

-0.524 

-1.349 

27 

0.150 

-0.496 

2S 

-0.460 

-1.350 

28 

0.000 

-0.495 

26 

-0.339 

-1.360 

29 

-0.180 

-0.496 

27 

-0.285 

-1.320 

30 

-0.346 

-0.480 

28 

-0.300 

-1.254 

35 

1.305 

-0.585 

29 

-0. 360 

-1  . 140 

36 

1.306 

-0.705 

30 

-0.383 

-1.071 

37 

1 .200 

-0.676 

31 

-0.376 

-1.007 

38 

1.050 

-0.676 

32 

-0.346 

-0.H84 

39 

0.900 

-0.676 

33 

-0.315 

-0.750 

40 

0.760 

-0.67*5 

34 

-0.299 

-0  • 70*> 

41 

0.600 

-0.676 

1 

-0.106 

0.000 

42 

0.460 

-0.675 

2 

-0.203 

0.000 

43 

0. 300 

-0.675 

3 

-0.300 

0.000 

44 

0.160 

-0.675 

4 

-0.323 

-0.060 

46 

0.076 

-0.600 

S 

-0.36S 

-0.180 

46 

-0.075 

-0.600 

6 

-0.411 

-0.287 

47 

-0.225 

-0.600 

7 

-0.461 

-0.405 

71 

1.320 

-0.825 

8 

-0.491 

-0.460 

72 

1.200 

-0.82-5 

9 

-0.460 

-0.496 

73 

1 .060 

-0.825 

10 

-0.390 

-0.670 

74 

0.900 

-0.82*5 

1 1 

-0.330 

-0.645 

75 

0.750 

-0.825 

12 

-0.285 

-0.698 

76 

0.600 

-0.825 

13 

-0.255 

-0.690 

77 

0.450 

-0.826 

14 

-0.158 

-0.690 

78 

0.300 

-0.826 

IS 

-0.060 

-0.690 

79 

0.150 

-0.825 

16 

-0.023 

-0.705 

80 

0.090 

-0.760 

17 

-0.015 

-0.750 

81 

1.320 

-0.960 

18 

-0.015 

-0.885 

82 

1 .200 

-0.960 

19 

-0.015 

-1.005 

83 

1 .050 

-0.960 

20 

-0.009 

-1.040 

84 

0.900 

-0.960 

21 

0.030 

-1.050 

85 

0.750 

-0.960 

22 

0.165 

-1.050 

86 

0.600 

-0.960 

23 

0.315 

-1.053 

87 

0.460 

-0.960 

24 

0.450 

-1.056 

88 

0.300 

-0.960 

25 

0.585 

-1.058 

89 

0.160 

-0 . 9 60 

26 

0.720 

-1.061 

90 

0.090 

-0.885 

27 

0.855 

-1.062 

32 

-0.600 

-0.670 

28 

0.990 

-1.064 

33 

-0.750 

-0.585 

29 

1.125 

-1.067 

34 

-0.900 

-0.600 

30 

1 .260 

-1.068 

49 

-0.460 

-0.676 

31 

1.365 

-1.070 

50 

-0.600 

-0.675 

32 

1.418 

-1 .050 

61 

-0.760 

-0.675 

33 

1.425 

-1.005 

52 

-0.930 

-0.690 

34 

1.425 

-0.870 

53 

-0.450 

-0.826 

3S 

1.425 

-0.735 

54 

-0.600 

-0.825 

36 

1.425 

-0.600 

55 

-0.750 

-0.825 

37 

1.421 

-0.525 

56 

-0.900 

-0.825 

38 

1.380 

-0.480 

57 

-1 .005 

-0.825 

39 

1.260 

-0.416 

58 

-0.480 

-0.975 

40 

1.140 

-0.344 

59 

-0.600 

-0.975 

4 1 

I .022 

-0.273 

60 

-0.750 

-0.975 

42 

0.902 

-0.207 

61 

-0.900 

-0.975 

*3 

0.797 

-0.144 

62 

-1.050 

-0.975 

44 

0.735 

-0.120 

63 

-0.480 

-1.125 

45 

0.660 

-0.117 

64 

-0.600 

-1.125 

46 

0.510 

-0.117 

65 

-0.750 

-1.125 

47 

0.360 

-0.119 

66 

-0.900 

-1.125 

48 

0.210 

-0.119 

67 

-1.050 

-1 .095 

49 

0.075 

-0.111 

68 

-0.450 

-1.245 

jo 

-0.036 

-0.066 

69 

-0.600 

-1.245 

1.440 

70 

-0.750 

-1.215 

f 0.20 

1.00 

2.35 

MMMM 


OUTPUT  DATA  - 1 34  - 


K > 2.35000  (l/FTI 


COMPLEX  VALUE  OF  DFDC  AT  THE  ENTRANCE  AND  AT  THE  COMMON  BOUNDARY  (1/FT) 

• l. 29793  -9.20573  -1.17960  -8.18918  -1.57852  -8.27092 

-1.09136  -5.86153  -1.16056  -6.60026  -1.71356  -9.60766 


ABSOLUTE  VALUE  OF  OFOC  AT  THE  ENTRANCE  (1/FTI 
9.29678  8.2737C 


ABSOLUTE  VALUE  OF  DFOC  AT  THE  CCMMON  BOUNDARY  (1/FTI 


8.42021 

5.94260 

6. 50463 

9.75925 

COMPLEX  VALUE 

OF 

THE  BOUNCARY 

F FUNCTION 

Q! 1 .11 

(BASIN  A) 

0.34629 

1.65656 

0. 34037 

1.39798 

0.36667 

1.64139 

0.63036 

2.11O07 

0.43959 

2.29322 

0.37951 

2.06259 

0.29358 

1.64680 

0.24941 

1.42348 

0.24769 

1.42556 

0.24661 

1.44330 

0.24957 

1.48275 

0.27716 

1.63593 

C. 36870 

2. 16C62 

0.44250 

2.61231 

0.49460 

2.97566 

0.48700 

2.9S361 

0.46333 

2.93952 

0.44983 

2.93968 

0.60990 

2.71710 

0.41472 

2.75151 

0.43526 

2.87409 

0. 37550 

2.50065 

0.26589 

1.83511 

0.  12728 

0.99263 

-0.03264 

0.01485 

-0. 19773 

-1.00370 

-0.35529 

-1.98789 

-0.48970 

-2.84220 

-0.58436 

-3.46495 

-0.63129 

-3.60009 

-0.62615 

-3.78956 

-0.60483 

-3.67437 

-0.66303 

-3.87777 

-0.63171 

-3.78176 

-0.60752 

-3.62347 

-0.58419 

-3.48298 

-C. 55838 

-3.32547 

-0.54655 

-3.23006 

-0.44466 

-2.58781 

-0. 31513 

-I.  79269 

-0.  17160 

-0.92932 

-0.03718 

-0.13393 

0.04134 

0.32326 

0.09129 

0.61347 

0. 18700 

1.1632  1 

0. 31682 

1.88821 

0.42339 

2.44775 

".479 22 

2.67833 

0.47367 

2.51123 

0.43217 

2.09534 

absolute  value 

OF 

THE  BOUNCARY  F FUNCTION 

(8ASIN  A ) 

1.49521 

1.63882 

1.68185 

2. 15351 

2.33497 

2.09721 

1.67277 

1.44516 

1.44692 

1.46422 

1.50361 

1.65924 

2. 19185 

2.64952 

3.01648 

3.03297 

2.97581 

2.97389 

2. 76785 

2.78259 

2.90686 

2.52869 

1.85427 

1.00076 

0.03585 

1.02299 

2.01939 

2.88408 

3.51388 

3.  85217 

3.84061 

3.72382 

3.93072 

3.83415 

3.67404 

3.53164 

3.37202 

3.27597 

2.62573 

1.82018 

0.94503 

0.13899 

0. 32587 

0.62022 

1.17816 

1.91460 

2.68610 

2.72087 

2.55551 

2.1 3*^45 

COMPLEX  VALUE 

OF 

THE  BCUNOARY 

F FUNCTION 

Ql l.ll 

(BASIN  ft  1 

0.24957 

1.68275 

0.24661 

1.44330 

0.24769 

1.42556 

0.26961 

1.62368 

0.18508 

1.07651 

0.08154 

0.69197 

-0.02073 

-0.C9970 

-0.09078 

-0.51546 

-0.11865 

-0.68622 

-0.12639 

-0.72387 

-0.15202 

-0.89069 

-0.18830 

-1.11269 

-0.22136 

-1.31929 

-0.23941 

-1.44184 

-0.22656 

-1.36291 

-0.22380 

-1.35925 

-0.24600 

-1.68099 

-0.24450 

-1.66753 

-0.26120 

-1.60105 

-0.24022 

-1.37715 

-0.23895 

-1.36033 

-0.23483 

-1.31758 

-0.23688 

-1.30970 

-0.23901 

-1.31001 

-0.26552 

-1.32932 

-0.23686 

-1.25921 

-0.23855 

-1.27718 

-0.23120 

-1.26359 

-0.19272 

-1.03660 

-0.16719 

-0.78037 

-0.07073 

-0.36280 

0.05222 

0.36133 

0.16735 

1.01662 

0. 16802 

1.01165 

ABSOLUTE  VALUE 

OF  THE  BOUNOARY  F FUNCTION 

(BASIN  SI 

1.50  361 

1.66622 

1.66692 

1.66516 

1.09230 

0.69868 

0.10183 

0.52339 

0.69640 

0. 73668 

0.90357 

1.12851 

1.33773 

1.66158 

1.38129 

1.37755 

1.50128 

1.66803 

1.62166 

1.39796 

1.38115 

1.33836 

1.33096 

1.33163 

1.  35180 

1.28092 

1.29926 

1.26689 

1.05636 

0.79613 

0.35002 

0.36508 

1.03030 

1.02551 

- 135  - 

HARBOR  OPENING  IFT.I«  0.200  OEPTH  (FT. I*  1.000 
K - 2.35000  ( I/FTI 

«»■  3.39500 

PERIOD  T « 0.72890  ISEC.I 


F2MAX«  3.9C599 


MESH 

PX 

PY 

FCMPLX 

F2 

FRA. 

l 

-0.210 

-0.105 

0.42423 

2.04361 

2.08717 

0.534 

2 

0.750 

-0.225 

0.02710 

0.24676 

0.24825 

0.064 

3 

0.600 

-0.225 

0.16902 

1.07050 

1.08376 

0.277 

4 

0.450 

-0.225 

0.30452 

1.83239 

1.85752 

0.476 

5 

0.300 

-0.225 

0.41714 

2.43353 

2.46902 

0.632 

6 

0.150 

-C.225 

0.48211 

2.72735 

2.76963 

0.709 

7 

0.0 

-0.225 

0. 49389 

2.69062 

2.73557 

C.  700 

8 

-0.150 

-0.225 

0.47472 

2.49220 

2.53701 

0.650 

9 

-0.270 

-0.225 

0.45247 

2.33993 

2.38328 

0.610 

10 

1.035 

-C.  405 

-0.32289 

-1.83039 

-1.85865 

-0.476 

11 

0.900 

-0.375 

-C. 17552 

-0.94066 

-0.95690 

-0.245 

12 

0.750 

-0.375 

-0.01827 

-0.00353 

-0.01661 

-0.005 

13 

0.60C 

-0.375 

0.13945 

0.91960 

0.93011 

0.238 

14 

0.430 

-0.375 

0.28466 

1. 74857 

1.77159 

0.454 

15 

0.300 

-0.375 

0.40116 

2.38653 

2.42000 

0.620 

16 

0.150 

-0.375 

0.47398 

2.74838 

2.78895 

0.714 

17 

O.C 

-0.375 

C. 49676 

2.80875 

2.85234 

0.730 

18 

-0. 150 

-0.375 

C. 47457 

2.62322 

2.66580 

0.682 

19 

-0.300 

-0.375 

0.42244 

2.30591 

2.34429 

0.600 

20 

1.200 

-0.525 

-0.50131 

-2.92889 

-2.97148 

-0.761 

21 

1.050 

-0.525 

-0.  38394 

-2. 19427 

-2.22760 

-0.570 

22 

0.900 

-0.525 

-0.23230 

-1.26624 

-1.28737 

-0.330 

23 

0.750 

-0.  52  5 

-0.06350 

-0.24966 

-0.25761 

-0.066 

24 

O.toOO 

-0.525 

0.  10664 

0.  75868 

0. 76614 

0.196 

25 

0.450 

-0.525 

0.26156 

1.65897 

1.67946 

0.430 

26 

0.300 

-0.52  5 

0.38474 

2.35326 

2.38450 

0.610 

27 

0.150 

-C.495 

0.46471 

2.75896 

2.79782 

0.716 

28 

0.0 

-0.495 

0.48740 

2.83180 

2.87344 

0.  736 

29 

-0.18C 

-0.495 

0. 43966 

2.49926 

2.53763 

0.650 

30 

-0.345 

-0.480 

0.34517 

1.94195 

1.97239 

0.505 

35 

1.305 

-0.585 

-0.57306 

-3.38578 

-3.43393 

-0.879 

36 

1.305 

-0.705 

-0.60161 

-3.56173 

-3.61218 

-0.925 

37 

1.200 

-0.675 

-0.54690 

-3.20432 

-3.25066 

-0.832 

38 

1.050 

-0.675 

-0.43143 

-2.47493 

-2.51225 

-0.643 

39 

0.900 

-0.675 

-0.2775C 

-1.52355 

-1.54862 

-0.396 

40 

0.750 

-0.675 

-0.10237 

-0.45771 

-0.46902 

-0.120 

41 

0.600 

-O.o75 

0.07556 

0.61037 

0.61503 

0.157 

42 

0.45C 

-0.675 

0.23778 

1.56969 

1.58759 

0.406 

43 

0.300 

-0.675 

0.36694 

2.31830 

2. 34716 

0.601 

44 

0.  150 

-0.675 

0.44810 

2.76858 

2.80461 

0.718 

45 

0.075 

-0.600 

0.47306 

2.84340 

2.88248 

0.738 

46 

-O.C  75 

-0.600 

0.46640 

2.74025 

2.77966 

0.712 

47 

-0.225 

-C.600 

0.37945 

2. 19967 

2.23215 

0.571 

71 

1.320 

-0.825 

-0.62905 

-3. 73627 

-3. 78869 

-0.970 

72 

1.200 

-0.825 

-0.57639 

-3.39036 

-3.43901 

-0.880 

73 

1.050 

-0.825 

-0.46387 

-2. 67032 

-2.71031 

-0.694 

74 

0.900 

-0.825 

-0.30952 

-1.70642 

-1.73426 

-0.444 

75 

0.  750 

-0.825 

-0. 13144 

-0.61157 

-0.62554 

-0.160 

76 

0.600 

-0.825 

C. 05078 

0.99499 

0.49758 

0.127 

77 

0.450 

-0.825 

0.21757 

1.49660 

1.51233 

0.387 

78 

0.300 

-0.825 

0.35135 

2.29143 

2. 31821 

0.594 

79 

0.150 

-0.825 

0.43873 

2.80463 

2.83874 

0.727 

80 

0.090 

-0.750 

0.46054 

2.88010 

2.91669 

0.747 

81 

1.320 

-0.960 

-0.64349 

-3.85262 

-3.90599 

-l.COO 

32 

1.20C 

-0.960 

-0.59489 

-3.51605 

-3.56602 

-0.913 

83 

1.050 

-0.960 

-0.48333 

-2.79269 

-2.83420 

-0.726 

84 

0.900 

-0.960 

-0.32810 

-1.81508 

-1.84449 

-0.472 

85 

0.750 

-0.960 

-0. 14808 

-0.70063 

-0. 71611 

-0.183 

86 

0.  600 

-0.960 

0.03528 

0.42068 

0.4221 6 

0.108 

87 

0.450 

-0.960 

0.20252 

1.43492 

1.44914 

0.371 

88 

0.300 

-0.960 

0.33697 

2.24635 

2.27149 

0.582 

89 

0.150 

-0.960 

0.42130 

2.75772 

2.78971 

0.714 

90 

0.090 

-0.885 

0.45028 

2.90394 

2.93864 

0.752 

32 

-0.600 

-0.570 

0.07017 

0.42361 

0.42938 

0.110 

33 

-0.750 

-0.585 

-0.04051 

-0.21789 

-0.22162 

-0.057 

34 

-0.900 

-0.600 

-0.11072 

-0.63643 

-0.64598 

-0.165 

49 

-0.450 

-0.675 

0.12541 

0. 75249 

0.76287 

0.195 

50 

-0.600 

-0.675 

0.02105 

0. 14474 

0. 14626 

0.037 

51 

-0.750 

-0.675 

-0.06691 

-0.37197 

-0.37794 

-0.097 

52 

-0.930 

-0.690 

-0.14112 

-0.81831 

-0.83039 

-0.213 

53 

-C.450 

-0.825 

C *00355 

0.06822 

0.068  31 

0.017 

54 

-0.600 

-0.825 

-0.06240 

-0.32821 

-0.33409 

-0.086 

55 

-0.750 

-0.825 

-0.12331 

-0.69709 

-0.70791 

-0.181 

56 

-0.900 

-0.825 

-0.16927 

-0.98180 

-0.99628 

-0.255 

57 

-1.005 

-0.825 

-0.19048 

-1.11804 

-1.13415 

-0.290 

58 

-0.480 

-0.975 

-0.11508 

-0.60975 

-0.62052 

-0.159 

59 

-0.600 

-0.975 

-0.14334 

-0.78854 

-0.80147 

-0.205 

60 

-0.750 

-C.975 

-0.17999 

-1.02290 

-1.03861 

-0.266 

61 

-0.900 

-0.975 

-0.21017 

-1.22209 

-1.24003 

-0.317 

62 

-1.050 

-C.975 

-0.22932 

-1.35796 

-1.37719 

-0.353 

63 

-0.480 

-1.125 

-0.19528 

-1.06415 

-1.08191 

-0.277 

64 

-0.600 

-1.125 

-0.20480 

-1.13769 

-1.15597 

-0.296 

65 

-0. 750 

-1.125 

-0.22114 

-1.25896 

-1.27823 

-0.327 

66 

-0.900 

-1.125 

-0.23631 

-1.37532 

-1.39548 

-0.357 

67 

-1.050 

-1.095 

-0.24143 

-1.43407 

-1.45425 

-0.372 

68 

-0.450 

-1.243 

-0.23324 

-1.27102 

-1.29224 

-0.331 

69 

-0.600 

-1.245 

-0.23346 

-1.29992 

-1.32072 

-0.338 

70 

-0.750 

-1.215 

-0.23077 

-1.31389 

-1.33400 

-0. 342 

